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Referat 

Die Dissipation von Energie von einem molekularen System in die Umgebung und die 
damit verbundene Zerstorung der Phasenkoharenz hat einen Einfluss auf mehrere physika- 
lische Prozesse wie Bewegung der Schwingungsmoden eines Molekiils, eines Ions in einer 
Falle oder einer Strahlungsfeldmode, sowic auf Exzitonen- und Elektroncntransfcr. Elek- 
tronentransfer spielt eine wichtige Rolle in vielen Bereichen der Physik und Chemie. 

In dieser Arbeit wird die Elektronentransferdynamik mit Bewegungsgleichungen fiir 
die reduzierte Dichtematrix beschrieben, deren Herleitung ausgehend von der Liouville- 
von Neumann-Gleichung iiber die Kumulanten-Entwicklung fiihrt. Durch Ankopplung an 
ein Warmebad werden dissipative Effekte beriicksichtigt. Zunachst wird diese Theorie 
auf Modellsysteme angewendet, um die verschiedenen Einfliisse der Umgebung auf De- 
population, Dephasierung und Dekoharenz besser zu verstehen. Dann wird die Dynamik 
von konkreten intramolekularen Transferreaktionen in realen Molekiilen berechnet und 
die Ergebnisse mit denen von Experimenten und anderer Theorien verglichen. Zu den 
untersuchten Systemen zahlen die Komplexe H2P— ZnP— Q und ZnPD— H2P. 
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Chapter 1 
Introduction 



The behavior of many quantum systems strongly depends on their interaction with the 
environment. The dissipative processes induced by interaction with the environment have 
a broad area of apphcations from isolated molecules to biomolecules. In order to achieve 
a realistic description of a molecular process, it is important to take the dissipation into 
account for systems like, e.g., vibrational levels in a big molecule, the quantized mode of 
an electromagnetic field, or a trapped ion. 

The interaction of a system with an environment also plays the main role in the mod- 
eling of electron transfer (ET) processes as it ensures their irreversibility. ET is a very 
important process in biology, chemistry, and physics |^, |, ||, || . It constitutes a land- 
mark example for intramolecular, condensed-phase, and biophysical dissipative dynamics. 
ET plays a significant role in nature in connection with conversion of energy. In the pho- 
tosynthetic reaction center, ET creates charge imbalance across the membrane, which 
drives the proton pumping mechanism to produce adenosine triphosphate. In chemical 
systems, surface ET between metals and oxygen is responsible for corrosion processes. In 
organic chemistry, mechanisms involving bond fracture or bond making often proceed by 
ET mechanism. In inorganic chemistry, mixed-valence systems are characterized by ET 
between linked metal sites. Finally, the nascent area of molecular electronics depends, 
first and foremost, on understanding and controlling the ET in specially designed chemical 
structures. That is exactly why the ET problem is the main topic of this work. 

Of special interest is the ET in configurations where a bridge between donor and ac- 
ceptor mediates the transfer. The primary step of the charge transfer in the bacterial 
photosynthetic reaction centers is of this type [Q, and a lot of work in this direction 
has been done after the structure of the protein-pigment complex of the photosynthetic 
reaction center of purple bacteria was clarified in 1984 [0. Many artificial systems, es- 
pecially self-organized porphyrin complexes, have been developed to model this bacterial 
photosynthetic reaction centers 0, B|, P]- The bridge-mediated ET reactions can occur via 



different mechanisms [T^, incoherent sequential transfer when the mediating 



bridge level is populated or coherent superexchange [|T^, 14 1 when the mediating bridge 
level is not populated but nevertheless is necessary for the transfer. In the case of the 
sequential transfer the influence of environment has to be taken into account. 

Apart from these aspects, one of the fundamental questions of quantum physics has 
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CHAPTER 1. INTRODUCTION 



attracted a lot of interest: why does the general principle of superposition work very well 
in microscopic physics but leads to paradox situations in macroscopic physics as for in- 



stance the Schrodinger cat paradox [15|. One possible explanation of the paradox and the 
non-observability of the macroscopic superposition is that systems are never completely 
isolated but interact with their environment [|1^. Interactions with the environment lead 



to continuous loss of coherence and drive the system from superposition into a classical 
statistical mixture. The question about the border between classical and quantum effects 
and systems, which model this problem, are also under considerations in this work. The 
interest in the decoherence problem is explained not only by its relation to the fundamen- 
tal question: "Where is the borderline between the macroscopic world of classical physics 
and microscopic phenomena ruled by quantum mechanics?", but also by the increasing 
significance of potential practical applications of quantum mechanics, such as quantum 



computation and cryptography |17, 18 



The rapid development of experimental techniques in the above-mentioned and other 
branches of physics and chemistry requires to describe, to model, and to analyse possible 
experiments by numerical and analytical calculations. The mathematical description of 
the influence of environment for all these examples has attracted a lot of interest but 
remains a quite complicate problem nowadays. The theory has been developed in recent 
years and a brief review of its progress is represented in the next chapter of this work. 
Despite of the intensive attention and investigations of this problem there is a necessity to 
consider basic model concepts in more detail in order to apply the mathematical techniques 
to real physical systems, which are studied experimentally in an appropriate way. 

In the present work we throw a glance on the known principles of the relaxation theory 
and we are mainly interested in the application of the relaxation theory to some concrete 
systems, such as a single vibrational mode modelled by a harmonic oscillator (HO), an 
artificial photosynthetic molecular aggregate, and a porphyrin triad, using simulations and 
numerical calculations as well as some analytical methods. The questions of the influence 
of environment on these systems are discussed in detail in this work. Our theoretical 
arsenal is based on the relaxation theory of dissipative processes containing calculations 
of coherent effects for electronic states and wave-packet dynamics, which provide the 
conceptual framework for the study of ET and decoherence. 

The basic concepts of the relaxation theory, based on the density matrix formalism, 
are reviewed in chapter ^ This technique will be used throughout this work. Chap- 
ter II deals with the question of the border between classical and quantum effects and 
reports on a study of the environmental influence on the time evolution of a coherent 
state or the superposition of two coherent states of a HO as a simple system display- 
ing the peculiarities of the transition from quantum to classical regime. Chapters ^ and 
^ concern the ET problem, namely the mathematical description of the ET in molecu- 
lar zinc-porphyrin-quinone complexes modeling artificial photosynthesis (chapter ^ and 
photoinduced processes in the porphyrin triad (chapter |^). Each chapter starts with an 
introduction and ends with a brief summary. The main achievements of the present work 
are summarized in the Conclusions. 



Chapter 2 

Reduced Density Matrix Method 



The goal of this chapter is to introduce the reader to the main mathematical tools for 
calculating the system dynamics induced by the interaction with an environment, which 
are used in all parts of this work. It should be noted that usually an environment is 
modelled by a heat bath, thus here and below we use environment and bath as syn- 
onyms for each other. The chapter starts with a review of historical developments of the 
theoretical models for dissipation processes and a brief consideration of various types of 
master equations, their characteristics, and different techniques used to describe these 
models (section |2.1| ). Using the Hamiltonian for system plus bath in the common form 
(section |2.2|) , the Green's matrix technique (section p.3|) , and the cumulant expansion 
method (section |2.4| ) we arrive at a differential form of the generalized master equation 
(GME) (secion 2.5), which is applied in the later chapters for the description of particular 
systems. Although some steps of this derivation are known in the literature we present 
them here in order to be rigorous and to reach completeness of notation in this work. 



2.1 Theoretical models for dissipation 



In the first quantum consideration of the system "atom + field" Landau [T^ has intro- 
duced an analog of the density matrix (DM), its averaging over the field states and its 
equation of motion. The rigorous introduction of the statistical operator and its equation 
of motion has been done by von Neumann in the early 30s |^. The first derivation 
of the master equations based on the Liouville equation for the system "atom -|- field 
-|- environment" were closely connected to the radio frequency range and to the satu- 
ration of signals in nuclear magnetic resonance 



21 



4[. For those systems the 
projector operator technique has been used in order to derive an exact integro-differential 
master equation |Q. One started to describe optical processes with the GME a bit 
later [||] b ecause in the 50s and earlier 60s there were still no experimental hints to the 



non-Markovian nature of the relaxation p6[ for optical processes. Nowadays a quantum 
dissipation theory has been much sought after as the goal of five communities: quan- 
^, |29|, ^ ^ , condensed matter physicists [^, | 



tum optics pi 



35| , mathematical 



physicists [l^l, astrophysicists gT]], and condensed phase chemical physicists 
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CHAPTER 2. REDUCED DENSITY MATRIX METHOD 



Theories of quantum dissipation can be divided into three main classes. The first 
class begins with a full system + bath Hamiltonian and then projects the dynamics onto 
a reduced subspace. Notable examples of this approach are the path integral approach 
of Feynman and Vernon Redfield theory ||2^, and the projection operator technique 
of Nakajima and Zwanzig ||^. As also mentioned by Pollard and Friesner |^ this 
theory can be divided into two subclasses: quantum and classical bath. At the opposite 
extreme the second class begins with linear equations of motion for the reduced density 
matrix (RDM) and then deduces a form for the equations of motion compatible with rele- 
vant characteristics of the relaxation theory. Examples of this approach are the semigroup 
approach of Lindblad pS] and the Gaussian ansatz of Yan and Mukamel |H6[ . An inter- 



and the Gaussian ansatz of Yan and Mukamel ^ ^ 

mediate procedure is to describe the bath as exerting a fluctuating force on the system. 
This approach is often used in the laser physics community, for example by Agarwal |2^ , 
Louisell 



30 



Gardiner |^ and others. This class of theory was also used for the exciton 
transfer description [Q. Pollard and Friesner denote this class as "stochastic bath 
theory" . 

In all mentioned theories the system is described with the RDM which evolves in time 
under master equations. It is known that the GME can be obtained by the following 
methods: (i) Decoupling of the relaxation perturbations and the DM of a full system 
25| , 0, ^ , (ii) averaging over the fast motions taking into account the hierarchy of char- 
acterisctic times ||5^, ^ 0, (iii) coupled multiparticle Bogolyubov equations (iv) 
Nakajima- Zwanzig projection operators technique (v) the diagrammatic tech- 

nique (vi) cumulant expansions (vii) the method of Green's functions for 

the full system, averaged over the realizations of the environment , and (viii) stochastic 
models [^, |6^ . 

The GMEs can be divided into two groups: integrodifferential GMEs (IGME), e.g., 
p3| , p8| , |60| , |66| derived using Nakajima-Zwanzig projection operators techniques or 
differential GME (DGME) [^, ^ often derived with the help of the cumulant ex- 
pansion technique 0, 0, |68|] . 

In most cases the IGMEs are transformed into the DGME at the outset, based on 

0' 



mm 



the assumption that the DM varies slowly over the bath correlation times [| 
59| , |66| . In linear optics the question of the relationship between IGMEs and DGMEs 
was discussed in |^, |7^. In the theory of relaxation it was shown (neglecting the effect 
of radiation on the relaxation) that IGME can also be constructed with the help of the 



cumulant expansion |7T 



It is generally assumed [^, |72[ that a fully satisfactory theory of quantum dissipation 
should have the following characteristics: (1) The RDM should remain positive semi def- 
inite for all time (i. e., no negative eigenvalues, which would in turn imply negative 
probabilities). Below we use the word "positivity" in order to point out this property. (2) 
The RDM should approach an appropriate equilibrium state at long times. (3) The RDM 
should satisfy the principle of translational invariance, if a coordinate and translation are 
defined. This condition requires that the frictional force be independent of the coordinate 
as is generally the case in the classical theory of Brownian motion. 

In the Markov approximation any GME can be reduced to one of the following types: 
Redfield (R), Agarwal (A), Caldeira-Leggett (C), Louisell-Lax (Lou), and Lindblad (Li). 



2.1. THEORETICAL MODELS FOR DISSIPATION 
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Now we simply list these types of master equations together with their characteristics. 



R. In the derivation of Redfield |^ one uses a quantum bath theory to model the 
environment. In the appropriate basis this master equation coincide with the equation 
of Agarwal's type. This equation is not of Lindblad form while in the Pollard and 
Friesner parametrisation the two types of master equations seem to be similar. The 
positivity of the DM which evolves under Redfield equation is violated for large system- 
bath coupling. At infinite time the DM reaches the thermal state of the bare system. The 
Redfield equation satisfies the requirement of the translational invariance. 

A. In the derivation of Agarwal's master equation |^ one models the environment with 
a stochastic bath. In the relevant basis this master equation coincides with the equation 
of Redfield type and does not coincides with Lindblad form. The positivity of the DM 
which evolves under Agarwal's equation is violated for large system-bath coupling. At 
infinite time the DM reaches the thermal state for the bare system. The Agarwal equation 
satisfies the requirement of translational invariance. System frequencies are modulated 
because of momentum and coordinate relax in different ways. The approach of Yan and 
Mukamel |^ ensures the same properties as theories by Redfield and Agarwal. 

C. In the derivation of the master equation of Caldeira-Leggett type ||53[ one models 
the environment with a quantum bath. One derives such a master equation with the path 
integral technique. A master equation of Caldeira-Leggett type is compatible with the 
equations of Redfield and Agarwal only in the high temperature limit. This equation is 
not of Lindblad form. The positivity of the DM is violated in this equation. The DM 
arrives at the equilibrium state only in the high temperature limit. The Caldeira-Legget 
master equation satisfies translational invariance requirement. 

Lou. In the Louisell-Lax approach one accounts for the environment with a 
quantum bath. It differs from Redfield and Agarwal equations by performing the rotating 
wave approximation (RWA) and agrees with Lindblad form. So the positivity of DM 
is maintained in this approach. At infinite time the DM reaches the thermal state of 
the bare system. In the equations of the Louisell-Lax type translational invariance is 
violated, because there is a coordinate dependent friction force. The system frequencies 
are constant in time. 



Li. The Lindblad master equation |3^ are constructed in a special form to conserve 
the positivity of the RDM. In the same basis this master equation coincides with the 
equation of Louisell-Lax type. At infinite time the DM reaches the thermal state, but 
the translational variance is violated. In order to preserve the translational invariance 
Lindblad 



38 



has included additional terms into the Hamiltonian and into the master 
equation. In that case the RDM at large times does not approach the equilibrium state 
of the bare system but some other state. This state is expected to be a projection of the 
equilibrium state of system plus bath onto the system subspace 

In our contribution we start with the demonstration of the derivation of the GME 
in differential form using the cumulant expansion technique. Afterwards, in chapter |] 
we arrive at a master equation of Agarwal type. In chapter ^ we use the GME in the 
Louisell-Lax form. The derived master equation with or without RWA is extensively used 
in all parts of this work. 
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CHAPTER 2. REDUCED DENSITY MATRIX METHOD 



2.2 Hamiltonian and density matrix 

In the common form the Hamiltonian "system + environment" can be written as 

H = H^ + H^ + H^^ (2.1) 

where H^ = I]-E^V^^ + hJ^'v^uVfiu represents a quantum system in the diabatic rep- 
resentation, E^ the energy of the diabatic state /i, H^ = J2 (j^t^^ + ^ bath of 
HOs, and H^^ = hJ2 {r^u + r+^,^ (v^^ + V^Jj the hnear interaction between them. Here 



= 1/^) (i^l is the transition operator of the system, v^j^y the couphng of the diabatic 
states /i and r^u = f'^tuh generahzed annihilation operator of the bath, 6f the 

annihilation operator of the bath mode ^ having the frequency u^, /C^^ the frequency- 
dependent interaction constant. Note, that rather often one factorizes the interaction 
constant on system and bath contributions. 



The DM of system plus bath p evolves under the von Neumann equation |^ p = 
—i/h [H, p]. The RDM a = Tr^p is obtained by tracing out the environmental degrees of 
freedom |7^. The coherent and dissipative dynamics of the system is described by the 
following equation 



-i/h 



H\ 



a 



H 



SE 



-i/h H^,a +Tr^(p) 



(2.2) 



where p denotes the DM p in the Heisenberg picture p = exp (i/fiH^tj p exp 
with respect to environment degrees of freedom. Then we substitute the unit operator of 
the form exp (—i/hH^t) exp (i/hH^t) before and after Tr^ (p) 



a 



-i/h 



H^,CT 



+ exp (-i/hHh) exp (i/hHh)Ti^ (p) exp (-i/hHh) exp (i/hHh). 

(2.3) 

Thus reduced density matrix equation of motion (RDMEM) takes the following form 



-i/h 



H^,a 



+ exp (^-i/hHh)TT^ (^p) exp {i/hHH), 



(2.4) 



where tilde denotes the interaction representation p(t) = exp (j/hH^t^pexp (^—i/hH^tj. 
Equation ( |2.4|) is in Schrodinger picture with respect to H^ and in Heisenberg picture 
with respect to H^; nevertheless it contains the interaction dynamics in the interaction 
picture. This representation is convenient because the relevant von Neumann equation 
for system plus environment p = —i/h H^^{t),p contains neither H^ nor H^. 



2.3 Green's matrix 



The commutator of an operator with the Hamiltonian can be represented symbolically as 
follows H^^{t),p = L(t)p where L{t) = L^^{t) is Liouville's operator in the interaction 



2.3. GREEN'S MATRIX 



15 



picture. Here we assume, that system and bath are disentangled at the initial moment of 
time 

p(0) = p(0) = p^(0)a(0). (2.5) 

In all calculations below we suppose that the initial states of the bath oscillators are 
thermalized pg(0) ~ exp{—hu^b'^b^/kBT). In accordance with |6^, ^] the evolution of the 
system is described by a Green's matrix 



p = DSE(t,0)p(0), 



(2.6) 



where D^^{t,0)p{0) = Texp 



-i/hJdTL{T) 



T stands for the time ordering operator. 



(2.7) 



The ansatz of the disentanglement Eq. (|2.5|) allows us to decouple the evolution of the 
system from the evolution of the environment on the basis of the reduced Green's matrix 
DSE(t,o) = Tr'^ [D^^{t,0)p^{0)] so that 

Tr^p(t) = Tr'^ [^^^(t, 0)p^(0)a(0)] = T>^''{t,0)a{0) 

The connection between a and p reads 

a = Tr^p = Tr^p = Tr^ [exp (^-i/hH^^pexp (i/hH^t 
= exp(^-i/hHH)Tr^pexp{i/hHh). 

Formula ( p.7| ) ensures the following property 

(D^^)-! (t, 0) Tr^p = a (0) . 

Substituting Eq. (|2.7| ) into Eq. ( p^) gives: 

a = -i/h [h^, a] + exp (-i/TiHh) t>^^{t, 0)(t (0) exp (i/hHh 



{2.t 



(2.9) 



(2.10) 



In the last equation the whole influence of the bath is included in D^^(t,0). In order 
to obtain a differential equation for a which is local in time we substitute the property 
Eq. (|]9D into Eq. (gJ^ so that 



a = ~i/h 



+ exp (-i/hHh) t>^^{t,^) 



D 



SEn-1/ 



(t,0)Tr^p exp U/kHH) . (2.11 



Factorizing the operator and the DM term one obtains 



a 



-i/h 

X 



+ exp {-i/hHH) ti^^it, 0)(DS^)-^(t, 0) exp {i/UHh 



exp (-i/hH^t) Tr^pexp (i/hHh 
Substituting Eq. (^ in Eq. ( pJ2D gives: 



(2.12) 



a 



-i/h 



+ exp (-i/hHh) ti^^it, 0)(D^^)-^(t, 0) exp (i/hHh) a. (2.13) 



So we have obtained a differential RDMEM instead of the integral Eq. 
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CHAPTER 2. REDUCED DENSITY MATRIX METHOD 



2.4 Cumulant expansion 



Equation ( p.l3| ) is found to be exact for the initially disentangled system and bath if 
the bath does not change in time. Such a precise description of the bath influence is 
possible with the path inetgral technique but found to be numerically expensive. Here 
we perform some approximations to consider the influence of the bath in leading order. 
Taking into account Eq. ( p.6| ) and Eq. ( |2.71 ) the reduced Green's matrix for Eq. ( |2.13| ) 
reads: 



DS^(t,0) = Tr'^|Texp J dr [-i/hL (r)] (0)| 



(2.14) 



Below we define a function F (A) as follows: 



Tr^ [F{\)] = (F(A)) = Tr^|TexpA ^ dr [-i/hL (r)] (0)| 



:2.15) 



To expand this expression in the cumulant form {F) = exp |x) A^i^nj we solve the 

following differential equation: ^F = —i/h\L(t)F. The solutions of this equation in 
zeroth, first, and second orders of perturbation theory are, respectively: 

t 

f(^) = 1 - xi/n J L{T)dT 



= 1 - \i/h J L (r) dr - Xi/h J L (r) dr -Xi/h J L (r') dr 



[2.16) 



The expression for D^^(t, 0) in second order of the perturbative and the cumulant expan- 
sions read, respectively 



(F(A,t)) ^ 1- AV^ / (L(r))rfr + 



-a 

n 



2 t T 



{L (r) drL {r')) dr' 







(F (A, t)) = exp 1^ X^kA ^l + Y. A"K„ + A"ir„) + . . . 
As Eq. (|2T7|) is equal to Eq. (12J8| ) so 



(2.17) 



(2.18) 



i 

Ki = -i/^J dr {L (r)) 



(2.19) 



1 x2 



t T 



^2 + 2 {KiY =[y) J J (r) drL (r')) dr' 





2.5. MASTER EQUATION 



17 



t T 







{L (r) L (r')) drdr' --\-i/h j {L (r)) dr 



[2.20) 



Now we shall use the well known fact that 



/ / (r) dr 



2jf{r)dTjf{T')dT'. This 





fact helps to rewrite the second cumulant Eq. ( 2.201) in the following form: 



= ( - ) jdrjdr'HL (r) L (r')) - {L (r)) (L (r'))} • 



(2.21) 







If we know Ki and K2 from Eq. ( |2.19| ) and Eq. ( p.21|) then the trace of F given by 
Eq. ( gjfSD yields: 

DS'^(t,0)= (F(A))|,^i=expK. (2.22) 

where K ^ Ki + K2. Equation ( p.22|) allows to express the kernel of the differential 
Eq. (|]T3D as DD^i = K, where 

d \ 
K = ^ (i^i + K,) = -z/n {L it)) + {-^/hf Jdr{{L {t) L (r)) - {L {t)) {L (r))} (2.23) 



Transforming back from Liouvillian to Hamiltonian form yields 

t 

(2.24) 



Ka 



;^^'^(t)),[(^^-(r)),a]]} 
It is easy to show that IH^^ (t) ) = since 



Ka = / rfr 



(2.25) 



This expression is found to be precise up to the second order cumulant expansion. Taking 
some approximation this equation can be transformed into either Agarwal-type master 
equation used for the analysis of the HO in chapter ^ or Louisell-Lax-type master equation 
used for the calculation of the ET ansfer dynamics in artificial photosynthetic molecular 
aggregates in chapter ^ Below we show some steps of derivation for the Louisell-Lax-type 
master equation. 



2.5 Master Equation 

Making the RWA in Eq. ( p.25[ ) and performing averaging over the bath degrees of freedom 
one gets bath correlation functions ^ \bx'{T)bx{t)) = dyx [n{ujx) + 1] exp [iuJx{t — r)], 
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where n{uj) = [exp {huj/ksT) — 1] ^ denotes the Bose- Einstein distribution. Back in the 
Schrodinger picture Eq. ( ^Hf ) reads 



X 

+ 



t/n[H'^,a]-l/h^T f' dr (2.26) 



where 



(2.27) 



are the correlation functions of the environment perturbations. The subscript n in 
Eq. ( |2.26| ) refers to the factor nx in Eq. ( ^^.271 ). To obtain Rn+i one replaces ra^ by + 1. 
The integral in Eq. ( p. 261) has different behavior on short and long time scales. On the time 
scale comparable to the bath correlation time the function R allows that non-resonant bath 
modes Ux ^ uo^y give a contribution to the system dynamics. Here we apply the Markov 
approximation, i.e., we restrict ourselves to the limit of long times and then the above men- 
tioned integral is an approximation of the delta function \imt^^ Rn{(jj,t — T)/nxdT = 
irY^x^x^i^ ~ ^x)- Furthermore, we replace the discrete set of bath modes with a 
continuous one. To do so one has to introduce the spectral density of bath modes 
J{uj) = TT J2x KlS{u — ux) and to replace the summation by an integration. Finally 
one obtains the following master equation 



a 



h 



[H'^,a]+La , 



with 



flU 



(2.28) 



(2.29) 



where the damping constant 

F^, = h-Mico^,) (2.30) 

depends on the coupling of the transition — > |z/) to the bath mode of the same 
frequency. Formally, the damping constant depends on the density of bath modes J at 
the transition frequency u^^. Equation ( p.29| ) belongs to Louisell-Lax type and maintains 
the Lindblad form. We apply this equation to the system of discrete levels to describe the 
ET process in chapter ^. A version of this equation without RWA is carefully investigated 
in chapter 3 in application to a single HO. 



Chapter 3 

First Application: Harmonic 
Oscillator 



In this chapter we develop and adopt the theoretical method, introduced in the previous 
chapter to the HO. This approach is well suited for systems with negligible electronic 
coupling between different diabatic states and a single reaction coordinate modelled by 
a HO. Although it is a quite simple system, the study of its dynamics allows to answer 
the question about the border between classical and quantum effects. This question deals 
with the superposition of coherent states of the HO. 

We begin this chapter with the introduction of the superpositional states and a brief 
review of decoherence problem (section In sections |3.2| and |3.3| we briefly rederive the 
methods of investigation of the HO coupled to a thermal bath. In section we discuss 
the behavior of the superpositional states either using the analytical method derived in 
section |3.3| or by numerical simulation. 

One of the goals of this contribution is to present a consistent analysis of the deco- 
herence on the basis of a DM approach starting from von Neumann's equation for the 
DM of the whole system, i.e. the microscopic quantum system and the "macroscopic" 
environment. 



3.1 Introduction to the decoherence problem 

There is a number of propositions how to create the superposition states in mesoscopic 
systems, or systems that have both macroscopic and microscopic features. A representa- 
tive example is the superposition of two coherent states of the HO 

\a,(f)) = N'^ (ja) + e^'^\-a)) (3.1) 

for a relatively large amplitude {a ~ 3^5). Here, |a) is a coherent state and = 
2 + 2 cos(/)exp ( — 2 |q;| ) is a normalization constant. These states have been observed 



recently for the intracavity microwave field [|76| and for motional states of a trapped ion 
||77| . Additionally, it has been predicted that superpositions of coherent states of molec- 
ular vibrations could be prepared by appropriately exciting a molecule with two short 
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laser pulses and the practical possibilities of realizing such an experiment have been 
discussed 



In this scheme the quantum interference would survive on a picosecond 
time scale, which is characteristic for molecular vibrations. 

From the theoretical point of view, quantum decoherence has been studied extensively 

Most efforts focused on the decoherence of the HO 



m 80, 81, 82, 83, m 



states due to the coupling to the heat bath, consisting of a large number of oscillators 
representing the environment. The energetic spectrum of the bath is usually taken to be 
broad and dense to provide the transfer of excitation energy from the system to the bath. 
The system is usually described on the basis of the master equation for the reduced density 
operator. There are few general approaches for this method. In most approaches listed 
in chapter ^ one adopts the Markov approximation for real calculations. It means that 
all details of the complex system-environment interaction are neglected and relaxation is 
described by the characteristic decay constants. The physical analysis of the system be- 
havior beyond the Markov approximation have been proposed by Zurek |T6[ : the coupling 
with the environment singles out a preferred set of states, called "the pointer basis" . Only 
vectors of this basis survive the quantum dynamics. The vectors of the pointer basis are 
the eigenvectors of operators, which commute with the (full) interaction Hamiltonian of 
the system. This basis depends on the form of the coupling. Very often this pointer basis 
consists of the eigenstates of the coordinate operator. The density operator describing 
the system evolves to diagonal form in the pointer basis, which is usually connected to 
the disappearance of quantum interference. The two approaches give different pictures of 
the same decoherence processes. 



3.2 Generalized master equation 



Let us consider a single molecule vibrational mode as a one-dimensional harmonic po- 
tential. For this case the common Hamiltonian Eq. (|2.1| ) contains the molecular system 
Hs = huj [a'^ a + 1 / 2) . The molecule interacts with a number of harmonic oscillators 
modeling the environment. In the interaction Hamiltonian 



HsE = nJ2K^{b+ + b^) (a + 



(3.2) 



a (a+) are annihilation (creation) operators of molecular vibrations with frequency u, 
6^ {bt) operators for the environmental vibrations having the frequencies u^. is the 



coupling between them. Performing the same steps of derivation as for Eq. ( p. 261) but 
without RWA we arrive at the non-Markovian master equation for the HO 



cr = —tuj 



a a, a 



La, 



(3.3) 



where the action of the relaxation operator L is defined by 



La 



(A„+i + A+) a, a+ + a] + [a+ + a,a (A„ + A+^^ 



(3.4) 
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Figure 3.1: Time-dependent relaxation coefficients 7n+i (solid line), (diamonds), 7n+i 
(boxes), 7„ (crosses) of the master Eq. ( |3.3|) , calculated for a bath containing 60 modes in 
the range [0,6w], with coupling function K{uj^) = VO.l. 



Here, the operators An and An+i are defined by the Unear combinations of the operators 
a and a"*" as 



An 
An+l 



7„ (t) a + 7„ (t) a+, 
7„+i {t) a + 7„+i (t) 



with the functions (see Fig. p.l|) 



7„+i(t) = / K+i(r)rfr = Kl{n^ + 1) 



7n(t) = f Rnir)dT = Y,Kln^- 



-i(co'g — C<j) 



7n+i(t) = / Rn+iir)e-^"'^dT = E^IK + 1. , 

JO ^ — ^(Co'^ + UJ 



%{t)= / Rn{T)e-^"'^dT = Y.Kln^- 



-i{uj^ + cj) 



(3.5) 

(3.6) 

(3.7) 
(3.8) 
(3.9) 



Note, that i?„ and Rn+i are the correlation functions of the environmental perturba- 
tions ( p.27[ ), where = [exp (^cu^/ZcsT) — 1]~^ denotes the number of quanta in the 
bath mode with the frequency u^. Subscripts n and n + 1 indicate the character of the 
temperature dependence, refering to the factors and + 1 in Eqs. (p.6|)-(pl9|). 
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The functions 7 correspond to the friction coefficient in the classical limit. The first 
two coefficients 7„ and 7„+i strongly depend on the coupling constant for frequencies 
a;^ ~ u; and on the number of quanta in the bath modes with the same frequencies, whilst 
the coefficients 7^ and 7n+i are very small for all frequencies. 

7„+i and 7„ describe the situation when an emission of a quantum from the system 
with rate 7n+i occurs more probably than an absorption of a quantum with rate 7^. The 
terms of the master equation associated with %+i and 7„ originate from the non-RWA 
terms ab^, a^ht of the Hamiltonian (|3.2|) and correspond to the reverse situation: an 



absorption 7„+i is more probable than an emission 7„. As shown in Fig. |3.1| the last two 
types of terms are essential only for the first stage of relaxation t <^ Tc, where Tc denotes 
the correlation time of environmental perturbations. 

The obtained master equation (|3.3| ) describes different stages of vibrational relaxation. 
The initial stage is defined by a period of time smaller than the correlation time r^. This 
time can roughly be estimated as Tc ~ 2tt/Auj from the width Auj of the perturbation 
spectrum K^. For such small times one can approximate the master equation 
in the following form 



cr 



— lUJ 



a] a, a + a 



+ 



a + a, 0" 



(a+ + a)]}, (3.10) 



where F = (2?t.^ + 1) is a real constant. As follows from Eq. ( |3.1CI| ), the pointer basis 

for this step of relaxation is defined by the eigenstates of the position operator Q ~ 



a. 



Another period of time, for which the form of the relaxation operator R according to 
Eq. ( |3.4| ) is universal, is the kinetic stage, where t ^ Tc and the Markov approximation 
becomes applicable. In this stage the master equation has the form 



(3.11) 



+7 {[[("- + l)a + 



na 



a + a, 0" 



\n + l)a + na^ | 



where 7 = TcK'^g is the decay rate of the vibrational amplitude. Here n = n^, K = 
and the density of the bath states g = giu)^) are evaluated at the frequency oj = oi 
the selected oscillator. It should be stressed that Eq. (|3.11| ) differs from the usual master 
equation for a damped HO for derivation of which the RWA is applied |76|, |79|, 0, 

82| , |83| , |8^ . Still Eq. ( p.ll| ) is only a particular case of the more general Eqs. ( p.3|) -( ^l9|) . 



To the best of our knowledge Agarwal was the first who derived this equation. The 
phase-sensitive relaxation leads to the new effect of classical squeezing and to a decrease 
of the effective HO frequency [^. Inbetween there is a time interval, where relaxation is 
specific and depends on the particular spectrum of K"^. 



3.3 Analytical solution for wave packet dynamics 



The solution of the equation of motion of the RDM can be conveniently found using 
the characteristic function formalism 



the differential operators ^ and ^ instead of a 



_d_ 

dx* 



89, 90, 91, 92 



This formalism enables us to use 
and a. From the set of normally 
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ordered F = Ti (^ae^'^^ e ^*°'^, abnormally ordered F = Ti (^ae ^*°-e^'^^^ and Wigner 

F = Ti ^(76^^*"+^"^^ characteristic functions we use here only the first. Multiplying 
both sides of Eq. (|3.3| ) with the factor / = exp {Xa~^) exp (— A*a) and taking the trace one 
can rearrange all terms into such a form that a and a"*" precede the appropriate exponent. 
For this operation we change the order of operators using the expression a exp (Aa+) = 
exp (Aa"*") (a + A) to make the transformation a"'"exp(Aa^) = ^exp(Aa"'"). After that 



every term can be represented by the normally ordered characteristic function 

F = Tr(/a), 



(3.12) 



upon which one of the differential operators acts. We obtain, e.g., Tr ([a+a, cr] /) 



d_ 



A^ F, Tr([aV,a]/) 



A* (1^ - A*) F, and Tr ([aa, a+] /) 



manipulations lead us finally to the complex-valued partial differential equation: 
P = _ [,^A* + /i (A + A*)] + M - /i* (A + A*)] 



-A*^F. Such 



(A + A*) (z/A* + z/*A) F, 



(3.13) 



where 



/x(t) = 7„+i(t)+7:(t) 



(3.14) 



are relaxation functions. An analogous method was used by Puri and Lawande |9^ , 
who also treated an HO coupled to the heat bath with the help of the normally ordered 
characteristic function. They obtained a general expression for time evolution of the char- 



acteristic function for an arbitrary initial state of the oscillator, see Eq. (12) in Ref. p2 |. 
This expression is valid for the non-Markovian regime but performed under RWA. 

We can solve Eq. ( |3.13D by using the integral representation for the characteristic 



function F (A, A*, 0) = Tr (^e'^^'^e "'^* j which formally allows us to describe the nondiagonal 

{a \P) is adopted. An initial characteristic 

"^-f^^' c{a, p) (3.15) 



DM. Below the notation / / dad(3c{a, (3) 
function 

F(A,A*,0) = j j dad(3e' 
will evolve in accordance with Eq. (|3.13|) as 



F(A, A*,t) = J J dadpexp 



c{a,p). 



(3.16) 



Taking the derivatives of the characteristic function ( p.l2|) one obtains the mean values 
of observables, in particular the mean number of quanta is given by 



a a 



d^F 



dXdX* 



dF dF 



A=A*=0 



(3.17) 



A=A*=0 



Here we restrict the cumulant expansion to the second order, i.e. m + n < 2 . For a 
wide class of initial states (coherent, thermal, squeezed, etc.) higher order cumulants 
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vanish and our approximation becomes exact. The cumulants could hold nondiagonal 
information, such as the DM, in relevant cases we stress it with the upper index (a) or 

The functions K^nif) in Eq. (|3.16| ) are given by the solutions of the sets of equations 

irfe) = (3.18) 



K20 = -u* + fiKn + 2{iuj- fi*)K2o, 
K02 = -u + fj.*Ku-2{iuj - fi)Ko2, 



(3.19) 



with the initial values 



i^n {t = 0) = K20 (t 



0) 
0) 
0) 



a, 



f3, 

K02 {t 



(3.20) 



0) = 0. 



3? and ^3 stands for real and imaginary part of a complex variable, respectively. For the 
special case P = a*, these initial conditions represent the coherent state with amplitude a. 
This solution can be used for the construction of wave packets in different representations. 
Here, we will discuss the coordinate representation, in particular the dependence of the 
probability density P on the vibrational coordinate Q and on time 



P{Q,t) 



2tt 



dAe-^^«x(A,t), 



(3.21) 



where 



x(A,t)=Tr 



a(a++a)^^^J = e-^'/^F {tX, -tX, t) 



(3.22) 



is a characteristic function for the position operator, which is nothing but the diagonal of 
the Wigner characteristic function [^. Evaluating the integral (|3.21| ) we finally obtain 



where 



and 



P{Q,t)= J J dad(3P'''''^^ {Q,t) c{a,(3), 



2^wjr) 



exp 



Q - Q("'^) (t) 



W{t) 



(3.23) 



(3.24) 



Q^-'^^ it) = it) + K^l> it) , 
V{t) = ^ + Kn (t) + K20 (t) + K02 (t) 



(3.25) 
(3.26) 
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For the case P = a* the function Q^°''°'*^ denotes the expectation values of the coordinate 
operator of the coherent state, V is the broadening of the Gaussian packet Eq. ( p.24|) . 
The distribution P can be used for the investigation of relaxation dynamics for any initial 
state of molecular vibration, but it is best suited for studying the evolution of the states 
prepared as a superposition of coherent states. Below, we will discuss the relaxation 
dynamics of initially coherent states and the superposition Eq. (|3.1|) of two coherent 
states. 



3.4 Coherent states and their superpositions 
3.4.1 Coherent states 

Here we investigate the relaxation dynamics using the master equation ( p.3|) . For the co- 
herently exited states |ao) the initial characteristic function is F (A, A*, 0) = exp («qA — ao\*) 
Thus, the initial values for k[q^ (t) and KqP (t) are and ao- In the first stage of re- 
laxation, when t <C Tc, the relaxation functions are /i [t) = 0, u {t) = Ft. Therefore, the 
solution of the system of Eqs. ( |3.18|) -( ^rT9D gives 



V{t) 



23fJ (aoe 



(3.27) 



Even in this early stage there is a small quadratic broadening of the wave packet pl"o'"oj 
{Q,t) without changing its mean amplitude. After the intermediate stage of relaxation 
the solution of the system goes into the Markovian stage of relaxation, where the master 
equation ( |3.11| ) works. For this stage the solution of Eqs. ( p.l8D -( P?TUD reads 

g(<^S,ao) 

V{t) 



where 



23ft [a^z it)] e"^* 
1 



+ 72 — ne 



-2-yt 



1 + 



7 



cos 2ojt) H — sin 2ujt 



(3.28) 



z [t) = cos out + {'-f /uj) sin ujt + i{u!/ uj) sin ujt, 

UJ = 



(3.29) 



n 



[exp{huj/kT) - 1] 



Equation ( |3.29| ) demonstrates the decrease of the effective harmonic oscillator frequency 
due to the phase-dependent interaction with the bath, well-known in classical mechanics, 
but absent in the majority of quantum considerations performed with RWA ^ 
94| , |95| . Neglecting the terms a'^bf , ab^ in Hamiltonian (|2.1|) , and repeating the derivation 



for the first-order cumulant moments one obtains the analogs of Eqs. (|3.18D 



k 



01 

(«) 

10 



{/3) 



- {lu + ^) K^l 



(3.30) 
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The solutions Kqi = exp {—iu — 7)t, KIq = exp {iu — 7)^ do not present any change 
of the HO frequency. So we conclude that this reduction of the frequency is due to the 
phase-sensitivity of the Hamiltonian ( |3.2| ). 




-6 -4 -2 2 4 6 8 
coordinate 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



variance 



Figure 3.2: Dynamics of the wave packet P{Q,t) (left plot), its variance l/(right plot, solid 
line), and mean number of quanta {a~^a) (right plot, dashed line) for 7 = 0.1a;, = Shuj, 
Qo = 4. 



The broadening of the wave packet V{t) is displayed in Fig. for the damping 
7 = O.lu chosen to present the prediction of Eq. ( |3.28| ) most clearly. An extremely 
rapid relaxation of the same order was assumed for the exciton motion simulation in a 
photosynthetic complex However, for the majority of real systems the damping rate 



is a few orders of magnitude lower. With methods of so-called quantum tomography 



even the rapid relaxation of the wave packet like in Fig. |3.2| can be observed experimentally. 
The increase of V{t) appears due to absorption of quanta from the heat bath and can be 
obtained already using RWA. The oscillation of the broadening which we directly derived 
in coordinate space is in accordance with the oscillations of second moments of the Green's 

Ref. 1^ assumes initially a squeezed state. Our prediction (|3.28|) 



function [27, 



99 



goes further, because the oscillation appears even for the usual coherent state as the 
initial one. This oscillation does not present quantum squeezing, because the width is 
never smaller than the ground state width. 

The oscillations with frequency 20 induced by the phase sensitivity of Hamiltonian ( p.l| ' 
also occur for the mean number of quanta (|3.17|) . For the Markovian stage of relaxation. 



a a 



n + e 



-2yt 



|«o| 



n 



■7 



— ] + 2|ao| — sin 2ct;t + n -:- I cos2Cjt 



7 



(3.31) 
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starts from the number of quanta for the coherent state |ao|^ and relaxes to the number 
of quanta in the resonant bath mode n, see Fig. |3.2| . Technically the oscillations are 
induced by the non-RWA relaxation functions 7^+1 and 7„, Eqs. (p.8|)-(p^, while Kohen, 



Marston, and Tannor treat the features of a non-RWA approach with the concept of 
a dissipation rate oscillating in time. 

In Fig. we have shown the dynamics of the mean number of quanta [a^a). Its 
distribution has been presented by Milburn and Walls | |100| | for a squeezed state of the 
HO. 



3.4.2 Superposition of two coherent states. Creation 



The methods derived in section are easy applicable to the evolution of the super- 
position of two coherent states. Superpositional states attract attention due to specific 
quantum effects they are involved in. Such effects of quantum nature which can be realized 
experimentally are discussed for quantum teleportation [ p.01| , p.02|| , quantum cryptography 
|TU3i , 



and quantum computation ||T^ |104| 



All these effects are possible as long as the 
states remain pure and keep their superpositional nature. 

Such states are created experimentally for the motional states of a trapped atom 
||77| , |105| | coupled to its hyperfine transition and for the microresonator mode interacting 
with a Rydberg atom ||106|| . The dynamics of such type of systems is successfully predicted 
by the so-called Jaynes-Cummings Model (JCM) ||107|| including a HO hojQa'^a coupled 
linearly to a two level system (TLS) with ^cutls via (in RWA) 



SB 



(3.32) 



where (T_ = |1) (2| denotes the TLS lowering operator, and allows the analytical solution 
1U7||. There are investigations considering the JCM coupling without RWA [ p.(J8| j. Another 

associated with 



109 



extension of JCM is the introduction of dissipation processes |^ 
spontaneous emission from a TLS and the loss of energy from the cavity through mirrors 



In the simplest generalization of the linear coupling the interaction is proportional 

Tn\. An intensity dependent interaction was considered by 
114|| has considered a multiquanta interaction. A particular case, 



to the number of quanta 
Buzek 



Gerry 



namely two-quanta transitions attracts interest up to nowadays ||1 14| , |115| , |116| , |117| , |118| . 

To illustrate the creation of superpositional cat-like states ( p.l| ) we have performed 
a calculation with the HO-TLS interaction Eq. ( |3.32| ). In our calculation the coupling 
strength g = 0.2c(Jtls is taken a few hundred times larger than in experiment ||106|| to 
make the wave packet dynamics more pronounced. 

We have calculated the time evolution of the DM pM^Nni^) of the TLS (M, de- 
noting, e. g., the Rydberg states of a Rb atom in [|106|| ) and a HO (m, n labeling the 
exitations of, e. g., the electromagnetic field mode in the microresonator which is tuned 
to cutls)- Initially the field is prepared in a coherent state, the atom in the excited state: 



f \ ^ 

PMmNnit = 0) = d2A/5jV2 exp -|a| 



(3.33) 
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The trace over the field degree of freedom yields the state of the atom 



a 



TLB 

MN 



TrnoP 



JCM 



MN 



Z^Pm 



JCM 

mNm 



and vice versa 



HO 



Pm' 



JCM 



rriMn ' 



(3.34) 



(3.35) 



m 



^mn — TrxLSP'^^^ - 

L J mn 

The field mode is characterized by the wave packet 

^"°(Q,i) = E^"°W-(Q)> (3-36) 
where ipmiQ) denotes the HO eigenf unctions, and the mean value of the coordinate is 



Q{t) 



P''''{Q,t)QdQ. 



The state of the atom is characterized by its population 



jTLS 



it) 



TLB 
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(t). 



(3.37) 



(3.38) 



The wave packet (|3.36|) shown in Fig. |3.3| (a) evolves from the coherent state with one 
peak to the superpositional state with two peaks (at uo^t ~ 27) or interference structure 
(at ujQt ~ 25). The mean value ( p.37| ) of the coordinate, see Fig. |3]^(b), refiects this 
transition from coherent to superpositional state: it oscillates initially and then relaxes to 
zero, because the wave packet of a superpositional state (|3.1| ) is symmetric with respect 
to Q = 0. The TLS population ( |3.38| ) in Fig. |^(c) also reduces to zero, because the 
interaction with each pair of HO modes n, n + 1 yields an oscillation of the TLS popu- 
lation (|3.38| ) with frequency given by the effective coupling g\/n + 1, so that the many 
HO levels initially populated, see Eq. ( |3.33|) , lead to destructive interference of these os- 



cillations. When the damping rate is small enough a revival occurs |^3|, |119|| . To prepare 
the superposition of coherent states we eliminate the revival by restricting the time of the 
atom- field interaction. Gerry and Knight ||120|| mention applying two 7r/2 pulses to the 
Rydberg atom before and after its interaction with a field mode. 

In this subsection we illustrated a possible method to create a superpositional state ( p.l| 
of a resonator field mode. The preparation of a motional state of a trapped ion in the 
state ( p.l|) is reported in Ref. while the excitation of the molecular vibration in this 

In the next subsection we 



state by two short laser pulses is discussed in Refs. 
calculate how the superpositional states evolve in time. 
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3.4.3 Superposition of two coherent states. Decoherence 

For the initial superposition of two coherent states ( p.l|) the normally ordered character- 
istic function Eq. ( |3.12| ) consists of four terms: 



F (A, A*, 0) = iV-2 + + e-2|°l' (e*<^F,.,_, + e-"t>F_^*,^ 



(3.39) 
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Figure 3.3: Preparation of superpositional state. The HO with ujq is coupled to a single TLS 
(f^TLS = l-Olwo), with coupling g = 0.2a;o- Initially the HO is prepared in a coherent state 
at Qq = 2.5, the TLS in the excited state. The time axis shows dimensionless time ujQt. (a) 
Evolution of the HO wave packet P{Q, t). (b) Mean value of the HO coordinate, (c) Population 
of the TLS. 

with Fa^p = e"'^"^''**. The first terms describe the mixture of two coherent states, |a) {a\ + 
|— a) (— a|. The last two terms correspond to a quantum interference, i.e., they reflect the 
coherent properties of the superposition. In accordance with Eqs. ( 3.18|) -( 2.21|) and the 



equation of motion (|3.11| ) of the RDM this initial state evolves as 

1 



P{Q,t) 
with the interference term 

Pint(Q,t) 



p("*>«) (Q, t) + p{-"*.-") (Q, t)] + Pi,t (Q, t) 



Ar2 



-2\a\' 



(3.40) 



(3.41) 



where P("'^) {Q,t) is given by Eq. QX^ ) with 

qK'") it) = -Q(-"*'-") (t) = 23? [az{t)] e-^\ 
g(-*'— ) (t) = -g(-°*'") (t) = 2i53 [az (t)] e-^\ 

V (t) and z (t) are defined by Eqs. ( ^.28| )-( [Og| ). Rewriting the real part of Eq. ( ^.411 ) we 
finally obtain 

I 2 



(3.42) 
(3.43) 



Pint(g,t) 



Ar2 



(3.44) 
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Figure illustrates, how the superpositional state ( |3.1|) evolves in time in accordance 




10 



time 



-2 2 
coordinate 

Figure 3.4: Time evolution of the superposition of coherent states for u; = 1, 7 = O.Olw, n = 0, 
a = 2, = f . 



with Eq. ( |3.40| ). It follows from Eq. (|3.44|) , that the interference term describing quantum 
coherence in the system is only significant when 



V{t) 



2 lal 



(3.45) 



This is true for 7t <^ 1, and moreover, when the two wave packets of the state come 
close together (i.e., at the moment when z (tj) ~ i). Expanding Eq. (|3.45| ) into a series at 
these points and taking into account 7^^ -C 1, it is easy to see that 



P,...~exp(-2|aP + »M^(l-2,y). 



(3.46) 



The decoherence is due to two reasons: The spreading V {ti) = 1/2 + AV {ti)n of the 
wave packet due to thermal excitations by the bath, and amplitude decoherence. The 
latter means, that even for the case n = and '^{az{ti)) = \a\ quantum interference 
disappears exponentially with the rate 



iJec^2|ar7 



(3.47) 



This result obtained by Zurek is the main reason why quantum interference is 
difficult to observe in the mesoscopic and macroscopic world. For example, a physical 



3.4. COHERENT STATES AND THEIR SUPERPOSITIONS 



31 



system with mass Ig in a superposition state with a separation of 1cm shows a ratio 
of relaxation and decoherence time scales of 10'^''. Even if our measuring device is able 
to reflect the quantum properties of the microsystem, nevertheless objectiflcation [ 121 ] 
occurs due to the coupling between the meter and the environment. The fundamental 
result of Eq. ( |3.47| ) is obtained no matter which approach is used, e.g. the RWA or Zurek's 
pointer basis approach or the self-consistent description of the present work. However, the 
time dependence of the superposition terms of the distribution of Eq. (|3.44| ) differs a little 



bit, which can be seen in Fig. 3.5. As quantum interference is more sensitive, we have used 



it for comparison of three different approaches to the present problem. Figure |3.5| shows 




time 



Figure 3.5: Time dependence of Pint {Q, t) for a; = 1, 7 = 0.25liJ, n = 0.4, a = 2, </) = 0, Qo = 0. 
The solid line represents the case Eq. (|3.28D , the diamonds represent RWA, the bullets represent 
earliest-time analysis Eq. ( p. 27 ). 



the difference between the time evolutions, which result from Eq. ( ^.271 ), from Eq. 
and from the corresponding result of the RWA approach. 

These differences arise due to the fact that during the relaxation there is no constant 
pointer basis for all steps of the evolution. As follows from Eqs. ( p.lOj ) and (|3.11| ), this 
basis changes from the position eigenstate basis in Eq. ( ^.10 ) to the more complicated 
basis in Eq. ( |3.11|) . 



3.4.4 Partial conservation of superposition 



As the interference term (|3.44|) of the wave packet decays inevitably, it is difficult to 
observe the effects mentioned in subsection p.4.2| . The decoherence determines, e. g., the 
main requirement to the potential elements of quantum computers: the decoherence time 
should be smaller than the computation time [1U4|. One should find a system maximally 
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isolated from an environment or increase the decoherence time artificially. One could 
prolong the coherent interval of the evolution by increasing the number of information 

fm\. 



transfer channels ||122|| , by feedback methods ||123|| , or by so-called passive methods 
Here we discuss in detail one more method of coherence preserving, namely organization 
of interaction with the environment |I6|, |125|, |126||. The idea of this method is to choose 



the system or the regime of system evolution which ensures a specific type of system-bath 
coupling, sometimes a nonlinear one. The linear coupling ( p. 21) to A = a'^ + a corresponds 
to the resonant exchange of quanta between system and bath. A coupling to the number 
of phonons used in description of a trapped ion |[L26|] and of exciton evolution in molecular 



crystals describes the bath-induced modulation of the system transition frequency. A 
coupling of the form A = {a+a){a — a) is discussed in Ref. ||125|| and Glauber's generalized 
annihilation operator ||127| , |128|| 



A 



aexp 



(3.48) 



in Ref. ||123|| . An exotic coupling operator a (a+a — (a'^a)) is discussed in Ref. ||126| . 



LU 
I- 

> 
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Figure 3.6: Schematic presentation of the two-quanta process. 



In Zurek's ansatz the eigenstates l'?/')^ of the coupling operator A (pointer states) 
are not perturbed by the interaction with the reservoir. For the linear coupling ( |3.2| ) and 
high temperature limit n ^ n + 1 in Eq. (|3.11|) such "pointer states" are represented 
by eigenstates of the coordinate operator [|16|, namely coherent states. In our approach 
they evolve in time. For the operator ( p. 4^ ) the eigen- and, respectively, pointer states 
coincide with the superpositional states |q;, vr) given by Eq. ( p.l| ). This coupling preserves 
the superposition from decoherence, but it is rather difficult to find an experimental 
system that provides such type of coupling to the environment. We consider a less exotic 
one, when the majority of the bath modes are not equal in frequency with the selected 
system, so that loss of amplitude and phase must be delayed. As a simplest example 
we take the quantum system surrounded by HOs with frequencies doubled to that of the 
system = 2uj as shown in Fig. 3^. Although this is still a resonance situation, it leads 
to unusual behavior compared to the usual case = u, discussed above. Describing it 
in RWA we rewrite the interaction (B.2|) 



H 



SE 



(3.49) 
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as discussed in Ref. [Q. Applying the evolution operator ( |2.7| ) we obtain again Eq. ( |3.3| ) 
for the RDM a of the selected system, but with some changes of the relaxational part 



La = T{n^ + 1) { 
+ rn^{[(a+)V, 



a 



+ 



} 



a^, cr 



(3.50) 



where F = irK^g^ is the decay rate of the vibrational amplitude. Here, the number of 
quanta in the bath mode n^, the coupling function K^, and the density of bath states 
are evaluated at = 2uj. An analogous equation was derived in Ref. ||129 
temperature limit was used in Ref. 



The zero 



In the basis of eigenstates \n) of the unperturbed oscillator the master eq. ( |3.50|) 
contains only linear combinations of terms as am,n = (^1 1^); o"m+2,n+25 and (Jm-2,n-2- 
It distinguishes even and odd initial states of the system. The odd state |1) cannot relax 
to the ground state |0), but the even state |2) can. 



o 
o 
o 




Figure 3.7: Influence of different baths. Evolution of the coherent state: Mean value of the 
coordinate, for Qq = —2.2. Diamonds: system coupled to the bath with lo^ = lo, ^ = O.lSw. 
Solid line: system coupled to the bath with uj^ = 2uj, T = 0.5uj. 
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The evolution of the system from different initial conditions was simulated numerically. 
The equations of motion of the DM elements are integrated using a fourth-order Runge- 
Kutta algorithm with stepsize control. To make the set of differential equations a finite 
one we restrict the number of levels by m, n < 20. 

We show the time dependence of the mean value of the coordinate in Fig. |3.7] . The 
mean value in the usual case uj^ = uj decreases with a constant rate. The same initial 
value of the system coupled to a bath with uj^ = 2uj shows a fast decrement in the first 
stage and almost no decrement afterwards. 

At temperature /c^T = 2huj/\n3), corresponding to n{2uj) = 0.5, we have simulated the 
evolution of the same superpositional states for uj^ = u and for lu^ = 2u, presented in 
Fig. ^]8|(a) and ^TsKb). For = u the quantum interference disappears already during the 
first period, while the amplitude decreases only slightly. The bath with = 2uj leaves 
the interference almost unchanged, although a fast decrease of the amplitude occurs. 
Therefore, it partially conserves the quantum superpositional state. 

3.5 Summary 

Starting from von Neumann's equation for a HO interacting with the environment mod- 
eled by a set of independent HOs we derived a non-Markovian master equation, which has 
been solved analytically. For two types of the bath with maximum of the spectral density 
near the system frequency and near the double of the system frequency and for two differ- 
ent initial states, namely a coherent state and a superposition of coherent states, the wave 
packet dynamics in coordinate representation have been analysed. It has been shown that 
wave packet dynamics demonstrates either "classical squeezing" and the decrease of the 
effective vibrational oscillator frequency due to the phase-dependent interaction with the 
bath, or a time-dependent relaxation rate, distinct for even and odd states, and partial 
conservation of quantum superposition due to the quadratic interaction with the bath. 
The decoherence also shows differences compared to the usual damping processes adopt- 
ing RWA and to the description using the pointer basis for decoherence processes. We 
conclude that there is no permanent pointer basis for the decay. There are two universal 
stages of relaxation: the coherence stage and the Markovian stage of relaxation, both 
having different pointer bases. We believe that the proposed method can be applied for 
other initial states and different couplings with the environment in real existing quan- 
tum systems, which is important in the light of recent achievements in single molecule 
spectroscopy, trapped ion states engineering, and quantum computation. 
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Figure 3.8: Influence of different baths, (a) Time evolution of the superposition of coherent 
states with initial separation 2Qo = 8 for o;^ = a; = 1, 7 = 0.005a;, n = 1.36. (b) Same as (a), 
but for = 2oj with V = 0.005a;, = 0.5. 



Chapter 4 

Electron Transfer via Bridges 



The present chapter deals with an important and quite difficult problem - the mathemat- 
ical modeling of ET. The purpose of our investigation is, at the one hand, to present a 
simple, analytically solvable model based on the RDM formalism [ p.3CI| , and to apply it 
to a porphyrin-quinone complex which is taken as a model system for the reaction centers 
in bacterial photosynthesis, and on the other hand, to compare this model with another 
one, which below we call the vibronic model. 

Before the description of these models and mathematical approaches some brief review 
of the ET problem in experimental and theoretical contexts is represented in section ^A . 
In section O we introduce the model of a molecular aggregate where only electronic states 
are taken into account because it is assumed that the vibrational relaxation is much faster 
than the ET. This model is referred to as the tight-binding (TB) model or model without 
vibrations below. The properties of an isolated aggregate are modeled in subsection 
4.2. 1| , as well as the static influence of the environment. The dynamical influence of bath 
fluctuations is discussed and modeled by a heat bath of HOs in subsection ^4.2.2| . The 
RDMEM describing the excited state dynamics of the porphyrin aggregate is described in 
subsection [4.2.3| and compared with an analogous equation of Haken, Strobl, and Reineker 
(HSR) p.31| , |132| , |133|| in appendix ^ In subsection [4.2.4] the system parameter 
dependence on the solvent dielectric constant is discussed for different models of solute- 
solvent interaction. In Subsection O system parameters are determined. The methods 
and results of the numerical and analytical solutions of the RDMEM are presented in 
subsection 

The dependencies of the transfer rate and final acceptor population on the system 
parameters are given for the numerical and analytical solutions in subsection [4.5. 1|. The 



analysis of the physical processes in the system is also performed there. In subsection 4.5.2 



we discuss the dependence of the transfer rate on the solvent dielectric constant for dif- 
ferent models of solute-solvent interaction and compare the calculated transfer rates with 
the experimentally measured ones. The advantages and disadvantages of the presented 
method in comparison with the method of Davis et al. [ffOl are analysed in subsection 



4.5.3. 



The vibronic model is described in section 4.6. In this case one pays attention to 



the fact that experiments in systems similar to the one discussed here show vibrational 
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coherence [ p.34| , |135|| . Therefore a vibrational substructure is introduced for each electronic 
level within a multi- level Redfield theory |p.36| , |137|| . The comparison of this model with 
the first one is done in section [4.7| . At the end of the chapter the achievements and possible 
extensions of this consideration are discussed. Unless otherwise stated SI units are used. 



4.1 Introduction to the electron transfer problem 

Long-range ET is a very actively studied area in chemistry, biology, and physics; both in 
biological and synthetic systems. Of special interest are systems with a bridging molecule 
between donor and acceptor. For example the primary step of charge separation in the 
bacterial photosynthesis takes place in such a system 



4.1.1 Mechanism of the electron transfer in bridge systems 



It is known that the electronic structure of the bridge component in donor-bridge-acceptor 
systems plays a critical role 



138, 13£ 



Change of a building block of the complex 
P, [L40| , |141|] or change of the environment [p.4CI| , |142|] can modify which mechanism is 
mainly at work: coherent superexchange or incoherent sequential transfer. Here the 
bridge energy stands for the energy of the state with electron localized on the bridge, 
not the locally excited electronic state of the bridge. When the bridge energy is much 
higher than the donor and acceptor energies, the bridge population is close to zero for all 
times and the bridge site just mediates the coupling between donor and acceptor. This 
mechanism is called superexchange and was originally proposed by Kramers [|143|| to de- 
scribe the exchange interaction between two paramagnetic atoms spatially separated by 
a nonmagnetic atom. In the case when donor and acceptor as well as bridge energies 
are closer than ~ k-^T or the levels are arranged in the form of a cascade, the bridge 
site is actually populated and the transfer is called sequential. The interplay between 
these two types of transfer has been investigated theoretically in various publications 
144| |145| , |146|| . Actually, there is still a discussion in the literature whether sequential 

or whether they are 



144 



transfer and superexchange are limiting cases of one process 
two processes which can coexist |^ . To clarify which mechanism is present in an artificial 
system one can systematically vary both energetics of donor and acceptor and electronic 
structure of the bridge. In experiments this is done by substituting parts of the complexes 



140| , |141| , |147|| or by changing the polarity of the solvent ||140|| . Also the geometry and size 



of the bridging block can be varied, and in this way the length of the subsystem through 



which the electron has to be transfered ||147| |148| , |149| , |150|| can be changed. 

Superexchange occurs due to coherent mixing of the three or more states of the system 
151| , |152|| . The transfer rate in this channel depends algebraically on the differences 
between the energy levels p, Q and decreases exponentially with increasing length of 
the bridge ||150| , |152|| . When incoherent effects such as dephasing dominate, the transfer 
is mainly sequential |Tl|, p.5CI|| , i. e., the levels are occupied mainly in sequential order 
0, |Tl], |12|, |140|| . The dependence on the differences between the energy levels is exponential 
H, M. An increase of the bridge length induces only a small reduction in the transfer 
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rate ||145| , |148| , |15(]| , p.52|| . This is why sequential transfer is the desired process in molecular 
wires 



150, 153 



4.1.2 Known mathematical theories 

In the case of coherent superexchange the dynamics is mainly Hamiltonian and can be 
described on the basis of the Schrodinger equation. The physically important results can 
be obtained by perturbation theory |l^, |154|| and, most successfully, by the semiclassical 
Marcus theory 



155 



The complete system dynamics can directly be extracted by nu- 
merical diagonalisation of the Hamiltonian ||150| , |156|| . In case of sequential transfer the 
environmental influence has to be taken into account. There are quite a few different ways 
how to include the influence of an environment modeled by a heat bath. The simplest 
phenomenological descriptions of the environmental influence are based on the Einstein 
coefficients or on the imaginary terms in the Hamiltonian ||157|, |158||, as well as on the 



Fokker-Planck or Langevin equations | p.57| , |158|| . The most accurate but also numerically 
most expensive way is the path integral method ||157|| . This has been applied to bridge- 
mediated ET especially in the case of bacterial photosynthesis | |159| ]. Bridge-mediated 
ET has also been investigated using Redfield theory [1^, |145| 

(e. g. [|5 



Liouville space |Tl 



by propagating a DM in 
In most of these methods 



and other methods (e. g. [ p.52| , p.56| , ^ 
vibrations are taken into account. 

The master equation which governs the DM evolution as well as the appropriate relax- 
ation coefficients can be derived from such basic information as system-environment cou- 
pling strength and spectral density of the environment |7^, |130| , |136| , |161| , |162|| . In the 
model without vibrations the relaxation is introduced in a way similar to Redfield theory 
but in site representation instead of eigenstate representation. A discussion of advantages 
and disadvantages of site versus eigenstate representation has been given elsewhere | 163| . 
The equations obtained are similar to those of Ref. |]TU| where relaxation is introduced in 
a phenomenological fashion but only a steady-state solution is found in contrast to the 
model used in this chapter. In addition, the present model is applied to a concrete system. 
A comparison of the ET time with the bath correlation time allows us to regard three time 
intervals of system dynamics: the interval of memory effects, the dynamical interval, and 
the kinetic, long-time interval |^3|]. In the framework of DM theory one can describe the 
ET dynamics in all three time intervals. However, often it is enough to find the solution 
in the kinetic interval for the explanation of experiments within the time resolution of 
most experimental setups, as has been done in Ref. |jlO|, |164|| . The master equation is an- 



alytically solvable only for simple models, for example ||158| , |165|| . Most investigations are 
based on the numerical solution of this equation |Tl], |136| , |137| , p.45| , |161| , |166|| . However, 
an estimation can be obtained within the steady-state approximation [l^, |167| . 

We should underline here that the RDMEM contains the coherent dynamics term 
and, therefore, is able to describe the SE mechanism. This fact has been recognized and 
successfully used to describe the reactions with dominance of superexchange in a various 
publications |Ty, 0, |145|| . Davis et al. |T^ have explicitly rederived the McConnel ana- 
lytical expression for superexchange using the RDMEM technique, based on the ability of 
the diabatic RDMEM to account for the coherent effects we apply the RDMEM formalism 
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to describe the ET in a real system where superexchange most probably takes place. 



4.1.3 Porphyrin-quinone complex as model for electron transfer 

The photoinduced ET in the supermolecule consisting of three sequentially connected 
molecular blocks, i. e., donor (D), bridge (B), and acceptor (A), is under consideration 
throughout this chapter. D is not able to transfer its charge directly to A because of their 
spatial separation. D and A can exchange their charges only through B. In the present 
investigation, the supermolecular system consists of free-base of tetraphenylporphyrin 
(H2P) as D, zinc-tetraphenylporphyrin (ZnP) as B, and p-benzoquinone as A [|140|| as 



shown in Fig. WA. 




Figure 4.1: Chemical structure of H2P — ZnP — Q. 



4.2 Model without vibrations 

In each of those molecular blocks shown in Fig. |4.1| we consider only two molecular orbitals, 
the lowest unoccupied molecular orbital (LUMO) and the highest occupied molecular or- 
bital (HOMO) |p.68|| . Each of the above-mentioned orbitals can be occupied by an electron 
or not, denoted by |1) or |0), respectively. This model allows to describe four states of 
the molecular block (e. g. D), the neutral ground state |1)homo|0)lumo (D), the neutral 
excited state |0)homo|1)lumo (D*),, the positively charged ionic state |0)homo|0)lumo 
(D+),, and the negatively charged ionic state |1)homo|1)lumo (D~),. Below a small ro- 
man index denotes the molecular orbital (m = - HOMO, m = 1 - LUMO), while a 
capital index denotes the molecular block (M = 1-D, M = 2-B, M = 3-A).A state 
of the supermolecule can be described as the direct product of the molecular block states. 

CMm = |l)Mm(0|Mm, CMm = |0)Mm(l|Afm, and flMm = cl^mCMm dcSCribc the CrcatioU, 

annihilation, and number of electrons in orbital Mm, respectively, while um = Z^m ^Mm 
gives the number of electrons in a molecular block. The number of particles in the whole 
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supermolecule is conserved J2m = const. Some of the electronic states of the molecular 



aggregate are shown in Fig. |472 . 




Figure 4.2: Schematic presentation of the energy levels in the H2P — ZnP — Q complex. The 
three states in the boxes play the main role in ET which can happen either sequentially or by 
a superexchange mechanism. Dashed lines refer to a sequential transfer, curved solid line to 
super exchange, dot-dashed to an energy transfer followed by charge transfer, dotted line optical 
excitation, and straight solid lines either fluorescence or non-radiative recombinations. 



Each of the electronic states has its own vibrational substructure. As a rule for the 
porphyrin containing systems the time of vibrational relaxation is found to be two orders 



of magnitude faster than the characteristic time of the ET ||169|| . Because of this we 
assume that only the vibrational ground states play a role in ET, and we do not include 
the vibrational structure. A comparison of the models with and without vibrational 



substructure will be given in section 4.6. 



4.2.1 System part of the Hamiltonian 

For the description the charge transfer and other dynamical processes in the system 
placed in a dissipative environment we use the common form of the Hamiltonian (|2.1|) 
where is the Hamiltonian of the supermolecule, the Hamiltonian of the dissipative 
bath, and H^^ describes their interaction. As mentioned in the introduction we are 
mainly interested in the kinetic limit of the excited state dynamics here. For this limit 
we assume that the relaxation of the solvent takes only a very short time compared 
to the timescale of interest for the system. Here the effect of the solvent is considered 
to be twofold. On the one hand the system states are shifted in energy. This static 
effect is state-specific and discussed below. On the other hand the system dynamics is 
perturbed by the solvent state fluctuations, which are independent of the system states. 
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The interaction Hamiltonian shall only reflect the dynamical influence of the fluctuations 
leading to dissipative processes as discussed in the next subsection. 

The static influence of the solvent is determined by the relaxed value of the solvent 
polarization and in general also includes the non-electrostatic contributions such as van- 
der-Waals attraction and short-range repulsion ||170| , |171|| . It is included into the system 
state energies and modeled as a function of the dielectric constant of the solvent. The 



static influence induces a change in the energy levels [|172|] , 



(4.1) 



where the energy of free and noninteracting blocks Hq corresponds to the energy of inde- 
pendent electrons in the fleld of the ionic nuclei. The term H^s denotes the state-selective 
electrostatic interaction within a molecular aggregate depending on the static dielectric 
constant eg of the solvent and V the inter-block hopping. It is assumed that the hopping 
V within the supermolecule is affected by the surroundings as discussed in subsection 



4.2.4. 



The energies of the independent electrons can be calculated by Hq = Y^Mm EMmnMrn, 
where denotes the energy of orbital Mm in the independent particle approximation 
1^, p.73|| . We introduce such a simplifled model to be able to calculate the energies of ionic 
molecular blocks e. g. D~. For each molecular block a fltting parameter Am is introduced 
which is used to reproduce the ground state-excited state transition e. g. D — > D*. Because 

the parameters Am are 
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these transitions change only a little for different solvents 
assumed to be solvent-independent. This is why we do not scale Hq. 



Emui one starts from fully ionized double bonds in each molecular block [173 



In order to determine 
calculates 



the one-particle states in the fleld of the ions in site representation and fllls each of these 
orbitals with two electrons starting from the lowest orbital; 



E 



Mm 



^ [Em- 2 Am cos 



2tcs\ 



(4.2) 



+ m Em - 2 Am cos 



27r 



M 



Here Nm and (Im are the total number of bonds and number of double bonds, respectively, 
in the porphyrin rings (M = 1, 2. For M = 3, i. e., the quinone (Q), see Sect. [4.3| .). The 
energy shift Em is chosen such that the neutral complex has zero energy. Am denotes 
the energy of hopping between two neighboring sites of the Mth molecular block. The 
function int() in Eq. (|4.2| ) denotes the integer part of a real number. In a similar way, 
by exciting, removing, or adding the last electron to the model system, one obtains the 
energy of the excited, oxidized, or reduced molecular block in the independent particle 
approximation. 

Below we apply Eqs. ( p.28|) -( ^l29D to the problem of evolution of a single charge- 
transfer exciton states in the system. In this case the number of states coincides with the 
number of sites in system: 

{Mm} /i. (4.3) 
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The next contribution to the system Hamiltonian is the inter-block hopping term 

(flM - 1)' + (flN - 1)'" . 
MN 

It includes the hopping operator between two LUMO states 

Vf^u = 4ncm, (4.4) 

as well as the corresponding intensities f^^^, i.e., the coherent coupling between different 
states of the system. We assume f 13 = because there is no direct connection between 
donor and acceptor. The scaling of v^y for different solvents is discussed in subsection 




Figure 4.3: Schematic view of the multiple cavities model 



The electrostatic interaction f^es scales like energies of a system of charges in a single 
or in multiple cavities surrounded by a medium with dielectric constant according to 
the classical reaction field theory ||174|| . Here we consider two models of scaling. In the 



first model each molecular block of the aggregate is in the individual cavity as shown in 
Fig. 4.3. For this case the electrostatic energy reads 

H,, = S''{e,)(H,i + H,on). (4.5) 



Here the function S^{es) describes the scaling of the electrostatic energy with the static 
dielectric constant eg of the solvent. The term 

H. = T.in,-l)^^ (4.6) 

takes the electron interaction into account while bringing an additional charge onto the 
block /i and thus describes the energy to create an isolated ion. This term depends on 
the characteristic radius of the molecular block. The interaction between the ions 

H:o. = T.Uf^.-m.-l)^^ (4.7) 

depends on the distance between the molecular blocks r^y. Both distances and r^y 
are also used in the Marcus theory | |155 |. The term + ifion reflects the interaction of 
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charges inside the aggregate which are compensated by the reaction field according to the 



Born formula ||175 | 



S 



H 



1 + 



1 



1^1 



(4i 



2es 2es 2 

In the second model sketched in Fig. 4.4, considering the aggregate as an single object 
placed in a cavity of constant radius one has to use the Onsager term ||175|| . This term 
is state selective, i.e., it gives a contribution only for the states with nonzero dipole 
moment, i.e., charge separation. Defining the static dipole moment operator as p = 



Y^ifi^j, — l){ny — l)f^,ye we obtain the Onsager term: Hes = 5 i^dip, where Hdip = — , 



S 



H 



2es 



(4.9) 




Figure 4.4: Schematic view of the single cavity model 

So, the physical meaning of both scalings is briefiy summarized as follows: Born- 
Marcus scaling corresponds to three cavities in the dielectric, each containing one molec- 
ular block of the aggregate. The energy scales with the Born formula ( |4.8| ). Onsager 
scaling Eq. ( [4.9|) refiects the naive idea that the whole supermolecule is placed in a single 
cavity. 



4.2.2 Microscopic motivation of the system-bath interaction and 
the thermal bath 



One can express the dynamic part of the system-bath interaction as 



(4.10) 



— * 

Here D^^{r) denotes the field of the electrostatic displacement at point r induced by the 
system transition dipole moment p^^, = p^iuiV^ + Vfj,u) ||158|| 
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(4.11) 



The field of the environmental polarization s denoted as P{f) = J2n ^{'^~^n)dn, where dn is 
the nth dipole of the environment and r„ its position. Only fiuctuations of the environment 
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polarization AP(r) influence the system dynamics. Averaged over the angular dependence 
the interaction reads 11172 



H 



SE 



E 
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(4.12) 



The dynamical influence of the solvent is described with the thermal bath model. 



The deviation A 



dr 



of dn from its mean value is determined by temperature induced 



dr 



But, solvent dipoles interact with each 



fluctuations. One could couple V^,^ to dn or A 

other even in the absence of a supermolecular aggregate. For unpolar solvents described 
by a set of HOs the diagonalisation of their interaction yields the bath of HOs with 
different frequencies ux and effective masses rrix- 

In the case of a polar solvent the dipoles are interacting rotators as, e. g. used to 
describe magnetic phenomena ||176| , |177|| . The elementary excitation of each frequency 
can again be characterised by an appropriate HO. So, we use generalized coordinates of 
solvent oscillators modes 



h 



2mxuJx 



(4.13) 



for polar as well as unpolar solvents. The occupation of the ith state of the Ath oscillator 
is defined by the equilibrium DM px^ij = exp (^—^^^6ij. 

All mutual orientations and distances of solvent molecules have equal probability. An 
average over all spatial configurations is performed. The interaction Hamiltonian ( |4.12|) 
is written in a form which is bilinear in system and bath operators: 
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SE 
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+ bx 



The coefficients 



2X^ 



d^7 



h 



S' 



SE/ 



(4.14) 



(4.15) 



depend on properties of the solvent, in particular, the frequencies ux- The precise de- 



termination of these coefficients needs special consideration. Expression ( [4.12| ) includes 
the dipole moment values corresponding to environmental modes and system transitions. 
The electric field of a dipole in medium is equivalent to the field of an imaginary dipole 
with a moment depending on the properties of the medium ||174|| . This influence of the 
medium is reflected here by the scaling function S^^. Explicit expressions for the solvent 
influence are still under discussion in the literature 1170, |171| . 



4.2.3 Reduced density matrix approach 

As usual the bath is given by HOs which describe the irradiative relaxation properties of 
the system. For the full description of the system one also should include photon modes 
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to describe for example the fluorescence from the LUMO to the HOMO in each molecular 



block transferring an excitation to the electro-magnetic field with rate [p.58|| . The rate 
of the radiative processes is small in comparison to other processes in system. That is why 
only the irradiative contribution is treated below. The treatment is similar to Redfield 



theory ||2^. For sake of notation and completeness we repeat the most important steps 
here for our model Hamiltonian. 

The irradiative contribution of the system-bath interaction corresponds to energy 
transfer to the solvent and spreading of energy over vibrational modes of the super- 
molecule. Applying the RWA to Eq. (|4.14| ) one gets 

^"'^ = EE^aPm^&^V + h.c, (4.16) 

where Kxp^,, denotes the interaction intensity between the bath mode 6a of frequency ujx 
and the quantum transition /xz/ between the LUMOs of molecules /z and v of frequency 
^iiu = I {Efi — Eu)- 

The dynamics of the system plus bath cannot be calculated due to the huge number 
of degrees of freedom. Therefore one uses the RDMEM technique. Below we assume that 
the coherent and dissipative dynamics can be represented by the independent terms of 
the RDMEM. This assumption is applicable if v^^ <C a;^,^. 

Here we use the RDMEM ( p.29|) derived in section p.5| . Here we treat the system 
Pfj,u and bath K\ contributions to the system-bath coupling separately. That is why we 
define the spectral density of bath modes as J{uj) = n J2x ^\^{^ ~ and the relaxation 
constant for Eq. ( |2.29D 

r^, = (4.17) 

depends on the coupling of the transition — > |z/) to the bath mode of the same 
frequency. Formally, the damping constant depends on the density of bath modes J at 
the transition frequency tu^jy and on the transition dipole moments between the system 
states p^u- 

For the sake of concrete calculations we write the RDMEM in matrix form. Substi- 
tuting the expressions for the exciton density a = J^o'fiu W) {^\, and operators Eq. ( [4.4| ) 

into the relaxation term Eq. ( p.29|) yields: 

+ r^A [n{uo^,x) + 1] + Tx^,n{uJx^)} (TnX- (4.18) 

A RDMEM of similar structure was used for the description of exciton transfer by Haken, 
Strobl, and Reineker in ||53|, [131| , |132| , |133| , |178|| . We present the comparison of these 



RDMEMs with Eq. ( |4T8D in appendix 0. 

For the description of a concrete system the introduced RDMEM ( p.28[ )-( p.29| ) is used 
in the matrix form Eq. ( 4.18|) with the use of an index simplification Eq. ( [4.3|) . For the sake 



of convenience of analytical and numerical calculations we replace the relaxation constant 
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T^i, and the population of the corresponding bath mode n{uj^y) with the intensity of 
dissipative transitions d^y = T^u\n{u!f^i^) \ between two states, as well as the corresponding 
dephasings intensity 



(4.19) 



With this one can express the RDMEM (|4l8| ) in the form 



'MM 



+ {La) 



a 



where 



{La] 



(4.20) 
(4.21) 

(4.22) 



From the manifold of system states we choose the ones which play the essential role in the 
electron transfer from the donor to acceptor Fig. |4.2| . They can be described in terms of 
single charge transfer exciton and correspond to the index simplification Eq. ( |4.3|) . The 
parameters controlling the transitions between the selected states are discussed in the 



next section and shown explicitly in Fig. E^S 
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Figure 4.5: Schematic presentation of the model without vibrational substructure. 



4.2.4 Scaling of the relaxation coefficients 

The relaxation coefficients F^,^ include the second power of the scaling function r^,^(e) = 
r°^S'^^ (e) because one constructs the relaxation term Eq. (|2.25|) with the second power 
of the interaction Hamiltonian. The physical meaning of H^^ is similar to the interaction 
of the system dipole with a surrounding media. That is why it is reasonable to use the 
relevant expressions from the relaxation theory of Onsager and Kirkwood ||174|| S^^ = 
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In the work of Mataga, Kaifu, and Koizumi ||179|| the interaction energy between 



2es + l' 

the system dipole and the media scales in leading order as 



2(er 



2es + l 



2eoo + 1 



(4.23) 



where eoo denotes the optical dielectric constant. In a recent paper of Georgievskii, Hsu, 
and Marcus ||170|| a connection J{uj) = —53a (a;) between J(a;) and the generalized suscep- 
tibility a{uj) is introduced and the explicit dependence a^u) on the dielectric parameters 
of the solvent is given. In application to the present work this means that the relaxation 
coefficient scales like r^^, = S^^ h''^J{uj^„)p'l^^. For example the approximation of spher- 



ical molecules gives a{uj) ~ ^ — ^. We approximate 
terms of scaling function it can be expressed as 



here, 



so 



i - ^. In 



SE 



1 
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(4.24) 



As an alternative possibility one can test the case of solvent-independent relaxation coef- 
ficient = const, S^^ = 1. 

The coherent coupling v^^, between two electronic states scales with eg and too, 
because a coherent transition in the system is accompanied by a transition of the envi- 
ronment state which is larger for the solvents with larger polarity. Here we involve the 
concept of reorganization energy from the model with vibrational substructure to account 
for this effect. For the reorganization energy we take the static part of the system-bath 
interaction calculated in frames of either individual or multiple cavity model. Here we 
take the static part H^^^ = S^^ (^H^i + -^ion) , note that the electronic part of H^^^ does 
not scale for the single cavity model. H^^'^ is associated with the so-called reorganization 
energy. The coherent couplings decrease with increase of H^^^. For the bath of relatively 
high frequency HOs (like the C 
as 



C stretching vibrations) this scaling can be taken [[18 



v% exp 



(4.25) 



where is the coupling of electronic states of the isolated molecule, tUvib the leading 
(mean) environment oscillator frequency. Unless otherwise stated c^vib = 1500 cm~^ is 
used. 



4.3 Model parameters 

The dynamics of the system is controlled by the following parameters: energies of system 
levels E^, coherent couplings f^i,, and dissipation intensities F^,^. The simplification is that 
we do not calculate these parameters, exept B state energy, rather take the corresponding 
experimental values. 

The absorption spectra of porphyrins ||140|| consist of a high frequency Soret band and 



a low frequency Q band. In the case of ZnP the Q band has two subbands, corresponding 
to pure electronic Q(0,0) and vibronic Q(1,0) transitions. In the free-base porphyrin 
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H2P the reduction of symmetry — ^ D2h due to the substitution of the central Zn ion 
by two inner hydrogens induces a sphtting of each subband into two, namely Q^{0,0), 
(5^(1, 0) and Q^(0, 0), 0). So the absorption spectra of ZnP and H2P consist of two 

and four bands respectively. In the emission spectra one sees only two bands for each 
molecule because of cascade intramolecular relaxations: pure electronic and vibronic one. 



Table 4.1: Low-energy bands of the porphyrin spectra in CH2CI2, from | 181 ] 
Absorption Emission 
Frequency, eV Width, eV Frequency, eV Width, eV 

~R2P i^So = 1-91 7oo = 0.06 z/o"i = 1-73 7oi = 0-05 
ZnP i/Qo = 2.13 7oo = 0.07 z/qi = 1.92 701 = 0.05 



Each of the above-mentioned spectra can be represented as a sum of Lorentzians with 
a good accuracy. The absorption spectrum of the complex H2P — ZnP — Q consists of the 
spectra of the individual H2P, ZnP, and Q compounds without essential changes. It means 
that intermolecular interactions in this case do not change the structure of electronic states 



of the subunits ||181| . We consider only the states corresponding to the lowest band of 
each spectrum to find the parameters of ET. The respective frequencies u and widths 7 are 
shown in Table T7L for CH2CI2 as solvent. The width of each band is formed by transitions 
to different sublevels of the electronic excited states of the aggregate and by the lifetime 
of these sublevels. In the present model consideration we do not have an excited level 
substructure, so that we cannot determine the linewidths correctly. Nevertheless we can 
consider the widths from Table |4.1| as upper limits for the relaxation constants. 

On the basis of the given spectral data and taking as reference energy -Edba = we 
determine Ed*ba = l/2(z/oo + ^01) = 1-82 eV and Edb*a = 2.03 eV (in CH2CI2) which 
allows us to extract the hopping energies Ad = 3.11 eV, A^ = 3.45 eV introduced in 
Eq. (|4.2|) . The excitation spectrum of the quinone lies far from the visible range, in the 
ultraviolet range. We do not have a precise spectral information about Q and, therefore, 
we do not calculate A a- We simply take the energy of the state with charge transfer to 



Q from reference [ p.4CI| |: i?D+BA- = 1-42 eV. | |182|| . Further Rempel and coauthors | ]14(j| 



estimate the electron coupling of the initially excited and the charged bridge states as 
(D*BA|if|D+B~A) = v^2 = 65 meV = 9.8 x 10^^ s~^ and they give the matrix element of 
two states with charge separation (D+B-A|if|D+BA-) = = 2.2 meV = 3.3 x 10^^ s"^ 
These values of the couplings are essentially lower than the energy differences between 
these corresponding system states 

hiUij > v°j. (4.26) 

This is the reason why it is useful to remain in site representation instead of eigenstate 
representation [|153|| . 

The relaxation constants are found with the help of the derived analytical formula to 
be r2i = r23 = 2.25 X 10^^ s~^. We discuss it in a more detailed form at the end of the 
next section. The typical radius of the porphyrin ring is about ri = r2 = 5.5 A, = 3.2 A 



141 1, while the distance r„j, between the blocks of the aggregate H2P — ZnP — Q reaches 



ri2 = 12.5 ± lA fng, r23 = 7 ± lA, = UA ± lA |TH3|, |nD|. 
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Table 4.2: Energy of the charged B state in different solvents and corresponding ET rates 
(for these calculations Born scaling (^]^) is used for the energies and Marcus scaling ( [4.24| ) 
for the relaxation constants). MTHF denotes 2-methyltetrahydrofuran, and CYCLO denotes 
cyclohexane. 



Solution 


CH2CI2 


MTHF 


CYCLO 


es, 181| 


9.08 


6.24 


2.02 


-Ed+b-a, eV 


3.12 


3.18 


3.59 


/cet, 10"^ s~^, num. 


33 


36 


0.46 


k^T, 10"^ s~^, an. 


33 


36 


0.46 


A;et, 10' s-\ exp. M 


23 ±5 


36 ±5 


+ 3 



The main parameter which controls ET in the triad is the energy of the state ii^D+B-A- 
This state has a big dipole moment because of its charge separation and is therefore 
strongly influenced by the solvent. Because of the special importance of this value we 
calculate it for the different solvents as a matrix element of the unperturbed system 
Hamiltonian Ho + Hes in Eq- (4.1). The calculated values of the bridge state D"'"B~A for 
some solutions are shown in Table ■ 



4.4 Results 



The time evolution of the ET in the supermolecule is described by Eqs. ( |4.2CI| )- (|4.21|) 



with the initial condition which corresponds to the excitation of the donor with a vr-pulse 
of appropriate frequency, i.e., the population of the donor is set to one. The system of 
equations was solved numerically and analytically. 

For the numerical simulation we express the system of Eqs. (|4.2CI|) - (|4.21|) in the form 



a = La, where a is a vector of dimension 3^ for the model with 3 system states and 
the super-operator L is a matrix of dimension 3^ x 3^. The time evolution of an element 
of the DM can be determined by a^uif) = cr^ti/(0) X^k^ W^mi^k? exp (A^gt), where A^^ and 
W^uK^ are the eigenvalues and eigenvectors of the super-operator L, respectively. These 
are obtained numerically. When fluorescence does not have to be taken into account, 
i.e., in the time interval t <^ G^~^ (cp. subsect. |4.2.3|) all states except |D*BA) (/x = 1), 



|D+B^A) (yU = 2), and |D+BA^) [ji = 3) remain essentially unoccupied, while those three 
take part in the intermolecular transport process. At later times t ~ G^"^ excitations 
decay into the ground state. So the physics of the ET processes in the molecular complex 
can be understood with the dynamics of the three above-mentioned states. The numerical 
simulation of the system dynamics with the parameters given in the previous section shows 
an exponential growth of the acceptor population. Such a behavior can be nicely fitted 
to a single exponential 

PsW = ^3(00) [1 - exp i-kETt)] , (4.27) 

where for the solvent MTHF A;et ^ 3.59 x 10® s'^ and P3(oo) ~ 0.9994. The population of 
the intermediate state /i = 2 which corresponds to charge localization on the bridge does 
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not reach a value of more than 0.005. This that means in this case the superexchange 
mechanism dominates over the sequential transfer mechanism. Besides it ensures the 
validity of characterizing the system dynamics with P3(oo) and 



A;ET = P3(oo)|^"[l-P3W]rft} • (4.28) 

The alternative analytical approach is performed in the kinetic limit 

t > l/min(7^,). (4.29) 

In Laplace space the inequality (|4.29| ) reads s <^ min(7^,^), where s denotes the Laplace 
variable. It is equivalent to replacing the factor l/iiuj^^ + '-f^^ + s) in the Laplace transform 
of Eqs. (|4.20|) - (|4.21|) with l/{iuj^^ + 7^^). This trick allows to substitute the expressions 
for non-diagonal elements of the DM Eqs. ( [4.21[ ) into Eqs. ( |4.20| ), so we do not use them 



explicitly anymore. After this elimination we describe the coherent transitions which 
occur with the participance of the non-diagonal elements by the following redefinition of 
the RDMEM (4.20)-(|;2|): 



cr 



MM 



9tMu(^tMtM + 9ufM(^UU- (4.30) 



In the given expression the transition coefficients g^u contain both, dissipative and coher- 
ent contributions ^ 

Qi^u = -I- v^,yVy^-f^y \h'^ (ujI^ + 7^^)] . (4.31) 

Now it is assumed that the bridge state is not occupied. This allows us to find the 
dynamics of the acceptor state in the form of Eq. ( |4.27| ), where 

, , 923(912 - 932) 

^ET = 923 H \ , (4.32) 

92\ + 923 

P,{oo) = ^''l'' (fcET)'^ (4.33) 

921 + 923 

The value of the dissipative coupling T^^ = h~^J{uj^,^)p'j^^, can be found by comparison 



of the experimentally determined ET rate and the derived analytical formula Eq. ( [4.321 ). 
To calculate J{uj^y) would require a microscopic model. We want to avoid a microscopic 
consideration and this is why we simply take the same T^i, for all transitions between 
excited states. The value of ET for II2P — ZnP — Q in MTHF is found by Rempel et al. 
140|| to be /cet = 3.6 ± .5 x 10^ s~^. If the bridging state has a rather high energy one 
can neglect thermally activated processes. V12 is negligible small with respect to ^23. In 
this case our result Eq. ( [4.32|) reads 



Tl UJ^i + i 21 + i 23 



Taking this formula, the relation between the dissipation intensities r2i = r23, and the 
experimental value for the transfer rate one obtains r2i = r23 — 2.25 x 10^^ s^^. The fit 
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of the numerical solution of the system of Eqs. ( [4.20| )-( P:.21| ) to the experimental value of 
the transfer rate in MTHF gives the same value. So the dissipative coupling constants 
are fixed for a specific solvent and for other solvents they are calculated with the scaling 
functions. With the method presented the ET in the supermolecule was found to occur 
with dominance of the superexchange mechanism with rates 4.6 x 10^ s~^ and 3.3 x 10^ s~^ 
for CYCLO and CH2CI2, respectively. 



4.5 Discussion 

4.5.1 Sequential versus superexchange 

Here we discuss how the transfer mechanism depends on the change of parameters. 
Namely which parameter have to be changed in order to alter not only the transfer rate 




11 12 13 14 




Figure 4.6: Dependence of the transfer rate (a) and the final population of acceptor state 
(b) on the coherent couplings vu (triangles and dashed line, V23 = ^23 = 2.2 meV), V23 (dots 
and solid line V12 = = 65 meV). Symbols correspond to numerical solution of system of 
Eqs. (|4.20D -(^l2T|). Lines correspond to the analytical result Eqs. (^I32|)-( ^l33|) . 



quantitatively, but the dominant mechanism of transfer and the qualitative behavior of 
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the system. In order to answer this question we calculate the system dynamic by varying 
one parameter at a time, while all other parameters are kept unchanged. The dependen- 
cies of transfer rate /cet and final population P3 on such parameters as coherent couplings 
V12, V23 and dissipation intensities r2i, are shown in the Figs. 4.6-4.7. The change of 
each parameter influences the transfer in a different way. 

In particular, /cet depends quadratically on the coherent coupling V12 from 10^^ s~^ 
to 10^^ s~^ in Fig. 4.6. Below it saturates at the lower bound /cet oc 3 x 10^ s~^. This 
corresponds to a crossover of the transfer mechanism from superexchange to sequential 
transfer. But, due to the big energy difference between donor and bridging state the 
efficiency of this sequential transfer is extremely low. This is displayed by P3(oo) ~ 0. 
In the region V12 ~ V23 both mechanisms contribute to the transfer rate. The transfer 
rate depends on V23 in a similar way. The decrease of final population in this region 
corresponds to coherent back transfer. At rather high values of V12, V23 — 10^^ s~^ the 
relation ( [4. 261 ) is no more valid. For this regime one has to use eigenstate instead of site 



representation because the wavefunctions are no more localized ||153|| . This variation of 
the coherent coupling can be performed experimentally by exchanging building blocks of 
the supermolecule. 




Figure 4.7: Dependence of the transfer rate on the relaxation intensities (triangles 
and dashed line), (dots and solid line). The rest of the system parameters corresponds 
to H2P — ZnP — Q in MTHF. Symbols correspond to numerical solution of the system of 
Eqs. ( 4.20 )-( 4.21 ), and lines to the corresponding analytical result Eqs. (|4.32| )-( 4.33 ). 



The variation of the dissipation intensities r2i, r23 in the region near the experimental 
values shows similar behavior of fcET(r2i) and /cET(r23) (see Fig. Here we assume an 
independent variation of r2i and r23. Both, fcET(r2i) and fcET(r23) increase linear until 
the saturation value 7 x 10^ s~^ at F > 10^^ s" 
between the numerical and analytical values. 



^ is reached. There is qualitative agreement 
In both dependencies one observes satura- 



tion at 10 s < F < 10 s . In Eq. ( 4.28 ) infinite time is approximated by 10 s 



It is found that approximating this time with larger values we do not obtain essential 
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changes in the results. But from the formal point of view we note that one cannot obtain 
transfer rates lower than 10^ s~^. Here we note that taking larger values for the infinite 
time does not change the transfer rates with the good precision. 

The physical meaning of the transfer rate dependence on the dissipation intensities 
seems to be transparent. At small values of F a part of the population coherently oscillates 
back and forth between the states. The increase of the dephasing Eq. ( [4.19| ) quenches 
the coherence and makes the transfer irreversible. So transfer becomes faster up to a 
maximal value. For the whole range of F, depopulations ^21? '^23 and thermally activated 
transitions di2, d^2 always remain smaller than the coherent couplings, therefore they do 
not play an essential role. 

Next, the similarity of the dependencies on F21 and F23 will be discussed on the basis 
of Eq. ( [4.32| ). Firstly, in the limit ksT /huj^i, the terms corresponding to thermally 
activated processes u^j^y < vanish and so \n{ujfj_y)\ = 0, while depopulations u^i, > 
remain constant \n{u^y)\ = 1. Secondly the condition u^^, ^ -^^y allows to neglect 7^^, in 
comparison with cj^^. With the above-mentioned simplifications the analytical expression 
for the transfer rate Eq. ( |4.32| ) becomes 

^ET ^ r2iF23 (F21 + F23)"' (ylJul^ + vlJul^ . (4.35) 

This equation is symmetric with respect to the relaxation intensities F21 and F23. This is 
why the transfer rate depends on each of them in the same way. 

The most crucial change of the transfer dynamics can be induced by changing the 
energies of the system levels. As was mentioned in section ^]3| these are the only pa- 
rameters which can be varied continuously by use of composed solvents. To the largest 
extent the mechanism of transfer depends on the energy of the bridging state |D"'"B~A). 
The results of the corresponding calculations are presented in Fig. 4.8. In different re- 
gions one observes different types of dynamics. For large energy of the bridging state 
-E'd+b-a ^ -E'd*ba and, respectively, E21 = -Ed+b-a ~ -Sd*ba ^ the numerical and ana- 
lytical results do not differ from each other. The transfer occurs with the superexchange 
mechanism. The transfer rate reaches a maximal value of 10^^ s~^ for low energies of the 
bridge. 

While the bridge energy approaches the energy of the donor state closer than ther- 
mal energy and goes even lower than this, the sequential mechanism of transfer starts 
to contribute to the reaction process. The traditional scheme of sequential transfer is 
obtained when donor, bridge, and acceptor levels are arranged in the form of a cascade. 
In this region the analytical solution need not coincide with the numerical result because 
the used approximations are no more valid. For equal energies of bridging and acceptor 
states k-ET displays a small resonance peak at E21 = —.4 as it is seen in Fig. 4.8(a). In 
the extremal case when the energy of the bridging state is even lower than the energy 
of the acceptor state a transfer does not take place anymore because the population gets 
trapped at the bridge state. The finite transfer rate for 

E21 < E31 (4.36) 

does not mean the actual ET because P3(oo) 0. For the dynamic time interval 
t < 7~J a part of the population tunnels force and back to the acceptor state with 
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the rate ^et- The analytical expression Eq. ( [4.32| ) gives constant rate for the regime 
( [4.36| ), while the numerical solution of Eqs. ([4.20|) - ([4.21|) displays instability, because such 
coherent oscillations of population cannot be described by Eq. ( 4.27| ) and /cet cannot be 
fitted with Eq. ( [4.28|) . In Fig. 4.8 the regime Eq. ( [4.36|) is displayed two times: for small 
E21 while E31 is kept constant and for large E31 while E21 remains constant. 
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Figure 4.8: Dependence of the transfer rate (a) and the final acceptor occupancy (b) on the 
energy of the bridging state E21 = -Ed+b-a ~ Ed*ba (triangles and dashed line, E31 = —0.4 eV) 
and the acceptor state E31 = -Eq+ba- ~ -£'D*ba (dots and solid line, £'21 = 1.36 eV). Symbols 
and lines correspond to numerical and analytical solutions, respectively. V12 = 65 meV, V23 = 
2.2 meV, T21 = T23 = 2.25 x 10^^ g-i 

The energy dependence of the final population has a transparent physical meaning for 
the whole range of energy. A large value of the bridging level ensures the transition of 
the whole population to the acceptor state with charge separation |D"*"BA~) which has 
the lowest energy of the excited states. In the intermediate case, when the bridging state 
has the same energy as the acceptor state, final population spreads itself over these two 
states ^3(00) = .5. 
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4.5.2 Different solvents 

Lowering the bridging state even more one arrives at the situation, where the whole 
population remains on the bridge as the lowest state of the system (taken into account 
here) and does not move to the acceptor. The dependence of the transfer rate on the 
relative energy of the acceptor E^i = -Ed+ba- — -Ed*ba in Fig. 4.8 remains constant 
while the acceptor state energy lies below the bridge state energy. Increase of i?3i up to 
E21 = 1.36 eV gives the maximal value of rate /cet oc r2i. While the value of £"31 increases 
further the acceptor level becomes the highest in the system and therefore the population 
cannot remain on it. 





— Born-Marcus, an. 

— Onsager-Mataga, an. 
X Experiment 
A Born-Marcus, num. 
♦ Onsager-Mataga, num 



Figure 4.9: The transfer rate k^T versus static dielectric constant. The energy of bridging 
and acceptor states scales in accordance with Born expression ( |4.8| ) (triangles and solid line), 
Onsager expression ( |4.9| ) (diamonds and dashed line) . Coherent couplings and dissipation scales 
in accordance with Mataga's expression ( [4.230 (diamonds and dashed line) , Georgievskii-Marcus 
expression ( |4.24|) (triangles and solid line). Symbols represent numerical solution of system of 
Eqs. ( 4.20| )- (|42l| ). Lines represents analytical result Eq. ( 4.32| ). Solid crosses with error bars 



correspond to experimental values [ 140 | 



For the application of the results to various solvents and comparison with experiment 
one should use the scaling for energy, coherent coupling, and dissipation as discussed 
above. The combinations of the energy scaling mentioned in subsection |4.2.1| and relax- 
ation intensities scalings mentioned in subsection 4.2.4 are represented in Fig. 4.9. An 
increase in the static dielectric constant eg from 2 to 4 leads to an increase of the trans- 
fer rate, no matter which scaling is used. Further increase of eg induces saturation for 
the Onsager-Mataga scaling and even a small decrease of the transfer rate for the Born- 
Marcus scaling. In all those cases the solvent is approximated as the continuous medium. 
Thus, the transfer rate depends on the interplay of the two mechanisms. Within the used 
approximations an increase in the solvent polarizability and, hence, of its dielectric con- 
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stant lowers the energy of the bridging and acceptor states and increases the system-bath 
interaction and, hence, the relaxation coefficients. It induces a smooth rise of the transfer 
rate in the whole interval of eg for the Onsager-Mataga scaling. On the other hand the 
large values of dielectric constant lead to essentially different polarisational states of the 
environment for the aggregate states with different dipole moment. The difference of the 
environmental polarizations reduces the values of the coherent couplings, see Eq. ( [4.25| ). 
This effect is reflected for the Born-Marcus scaling in the small decrease of /cet for large 
values of eg. The values for the transfer rate, obtained with this scaling come closer to the 
experimental value k-^Ti^^^^'"^^) ■ This gives a hint that the model of individual cavities 
for each molecular block is closer to the reality than the model with a single cavity for 
the whole supermolecule. 

Below we consider Born scaling Eq. ( |4.8| ) for the system energies and Marcus scaling 
Eq. ( |4.24| ) for the dissipation parameter to compare the calculated transfer rates with the 
measured ones. For the solvents CYCLO, MTHF, and CH2CI2 one obtains the following 
relative energies of the bridging level E21 = 1.77 eV, 1.36 eV, and 1.30 eV, respectively, 
i.e., a decrease in the energy of the bridging state. 

The calculated transfer rate coincides with the experimental value | 140 | for H2P— 



ZnP— Q in CYCLO, see table |4.2| . For CH2CI2 the numerical transfer rate diverges 
from the experimentally determined one. The calculated numerical value is found to be 
approximately thirty percent faster. The following reasons could be responsible for this 
difference: (i) absence of vibrionic substructure of the electronic states in the present 
model; (ii) incorrect dependence of system state energies on the solvent properties; (iii) 
appearing of additional transfer channels not mentioned in the scheme shown in Fig. 4.2. 
Each of these possibilities requires some comments. 

ad (i): The incorporation of the vibrational substructure will result in a complication 
of the model with parameters such as the mass of each vibrational mode [0, |180|| . It 
will also give a more complicated transfer rate dependence on the energy of the electronic 
states and dielectric constant. In the model with vibrational substructure the interplay 
between the difference of the energies of each pair of electronic states and corresponding 
reorganization energy determines if the pair belongs to the normal, activationless, or 
inverted region according to Marcus ||155|| . In contrast to the present consideration the 



model with vibrational substructure should yield the maximal transfer rate for nonequal 
energies of electronic state, namely for the activationless case: the energy difference equals 
to the reorganization energy. For a detailed comparison of TB model and the model with 
vibrational substructure see the next section. 

ad (ii): In particular such solvent effect as solvation shell, see for example Ref. | |184|| , 
do need a molecular dynamics simulation. The total influence of solvent is, probably, 
reflected in solvent-induced energy shift between the spectroscopically observable states 
without charge separation -E'd*ba and -E'db*a ||140|| which is neglected in our consideration. 

ad (iii): Some states of our three-site system have not been included in the model 



schematically presented in Fig. [4.2| . Using the solvent with strong dielectric constant can 
bring high-lying states closer to the ones involved in the transfer. The state which might 
play a role is |D^B+A) because of its larger dipole moment. So it is strongly influenced 
by the solvent. The state |DBA*) which has such a high energy that one can excite it only 
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with ultraviolet radiation, most probably, does not play a role in the discussed transfer. 



4.5.3 Comparison with the steady-state solution 

In the work of Davis et al. |T^ the vibrational substructure of the electronic states is 
also not taken into account, the relaxation is incorporated phenomenologically. We use a 
similar approach in the present consideration derived within a Redfield-like theory. Also 
we consider the relevant processes of dephasing and depopulation between each pair of 
levels. In contrast Davis et al. apply relaxation only to selected levels. In their paper 



dephasing 7 occurs between excited levels, see Eq. (7) in Ref. |jTO[, while depopulation 
k takes place only for the transition from acceptor to ground state. The advantage of 
the approach of Davis et al. is the possibility to investigate the transfer rate dependence 
for more than one bridge state. This was not the goal of the present work but it can be 
extended into this direction. We are interested in the ET in a concrete molecular complex 
with realistic parameter values and realistic possibilities to modify those parameter. Our 
results as well as the results of Davis show that ET can occur as coherent (with the 
superexchange mechanism) or dissipative process (with the sequential mechanism). We 
have considered various ways of switching between these two mechanisms including the 
one suggested by Davis et al., i. e., the change of the relaxation intensities meaning change 
of the solvent. The numerical steady state method used by Davis et al. is an attractive one 
due to its simplicity, but unlike our method it is not able to give information about the 
time evolution of the system. On the other hand the method of Davis et al. does not allow 
to calculate the most widespread stationary characteristics that clarify the presence of ET 
process, namely the fluorescence of the donor and its quenching. We suppose this can 
be included into the approach by Davis et al. by introducing the depopulation coefficient 
describing the transfer from donor to ground state. 



4.6 Vibronic model 

In the model with vibrational substructure (VSS) the influence of the interaction with 
an environment onto ET and other dynamical processes is also described by the common 
Hamiltonial of the form Eq. (|2.1|) . 

The bath again is modeled by a distribution of HOs and characterized by its spectral 
density J{uj). Starting with a DM of the full system, the RDM of the relevant (sub)system 



is obtained by tracing out the bath degrees of freedom |^ . While doing so a second-order 
perturbation expansion in the system-bath coupling and the Markov approximation have 
been applied |^. So one arrives at the RDMEM of the already discussed structure. 

The bridge ET system H2P — ZnP — Q is modeled by three diabatic electronic po- 
tentials, corresponding to the states |1) = |D*BA), |2) = |D+B-A), and |3) = ID+BA") 
(see Fig. [4.10|) . Each of these electronic potentials has a vibrational substructure. The 
vibrational frequency is assumed to be 1500 cm~^ as a typical frequency within carbon 
structures. The potentials are displaced along a common reaction coordinate which rep- 
resents the solvent polarization [ |155| j. Following the reasoning of Marcus |155|| the free 
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Figure 4.10: Electronic potentials and parameters of the vibronic model. The donor surface 
H2P* — ZnP — Q is given by the solid line, the bridge 112?"^ — ZnP~ — Q by the dashed line, 
and the acceptor H2P'^ — ZnP — Q~ by the dotted line. 



energy differences AG^j^ corresponding to the ET from molecular block u to fi {u = 1, 
= 2, 3) are estimated to be [|||, |l40l 

AG^n = - - - + AG^^ie,) (4.37) 

47reoes r^^ 

with the term AG^jy(es) correcting for the fact that the redox energies E^ and El'^'^ are 
measured in the reference solvent with dielectric constant e^''^: 

The excitation energy of the donor H2P —>■ H2P* is denoted by E^^. denotes the radius 
of either donor (1), bridge (2), or acceptor (3) and r^^ the distance between two of them 
as presented in section 0[ 



Also sketched in Fig. |4.1CI| are the reorganization energies A^^, = A^^ + A^^. These 
consist of the internal reorganization energy AJ^^, which is estimated to be 0.3 eV ||140|| , 
and the solvent reorganization energy ||155|| 

-2/1 1 1 \ / 1 1 



Further parameters are the electronic couplings between the potentials. They are the 
same as described in section Ol The damping is described by the spectral density J{uj) 
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of the bath. This is only needed at the frequency of the vibrational transition and is 
determined J(a;vib)/^^vib = 0.372 by fitting the ET rate for the solvent methyltetrahydro- 
furan (MTHF). In the vibronic model the spectral density is taken as a constant with 
respect to eg. 

Next the calculation of the dynamics is sketched. Starting from the Liouville equation, 
performing the abovementioned approximations the equation of motion for the RDM pMAf 
can be obtained [|136| , |137|] 

Q^PMAf = ^i^M - Ej^)PMAr - i "^{vaticPmic - vicmPjcn} + Rmn ■ (4.40) 

The index M. combines the electronic quantum number /i and the vibrational quantum 
number M of the diabatic levels Em- vmn = V^yF^c^p, M.v, N) comprises Franck- 
Condon factors Fpc and the electronic matrix elements V^u- The third term RmN de- 
scribes the interaction between the relevant system and the heat bath. In principle one 
can eliminate the coherent terms containing V]cm diagonalising the Hamiltonian of the iso- 
lated system. In such approach one couples the environment transitions to the transitions 
between the eigenstates of the system. This is adiabatic approach, used in e.g. ||145|| . It is 
rather expensive numerically. We do not apply it here, so we use the diabatic approach. 
There are still a discussion in the literature wether the diabatic method is able to provide 
precise results ||185|| . Nevertheless one uses the diabatic approach rather often ||136| , |137| 



because it is less expensive numerically ||186|| . 

Equation ( |4.40| ) is solved numerically with the initial condition that only the donor 
state is occupied in the beginning. The population of the acceptor state 

^3(t) =EP3M3Af(t) (4.41) 

M 

and the ET rate given by Eq. ( [4.28| ) are calculated by tracing out the vibrational modes. 



4.7 Comparison of models with and without vibra- 
tions 

Here we compare the following models: (i) the model where only electronic states without 
vibrational substructure are taken into account (see Fig. |475| ) and relaxation processes 
take place between the electronic states and (ii) the model where the relaxation takes 
place between vibrational states within one electronic state potential surface introduced 
in the previous section. 

In this section the energies of the electronic states Em of model without vibrational 
substructure are chosen to be the ground states of the harmonic potentials of the vibronic 



model shown in Fig. [4.10| . So they vary with the dielectric constant. The electronic 
couplings are chosen to scale as Eq. ([4.25|) . Note, that in the model with vibrational 
substructure it corresponds to the scaling provided by the Franck-Condon overlap elements 
between the vibrational ground states of each pair of electronic surfaces 

\\ 



VMM = V^yFpc (/^, 0, z^, 0) = V^, exp . (4.42) 
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Here we have used the expression for the Franck- Condon factor for the parabolas with 
equal effective mass m, and equal curvature FpcifJ^, 0,1^,0) = exp ( — '"^^"^ '^ifiu) def- 
inition of reorganization energy A^,y = mujl^^^^q^^^. q^u stands for the difference of the 
reaction coordinate equilibrium points of these parabolas. In the vibronic model not only 
the free energy differences AG but also the reorganization energies A scale with the di- 
electric constant e^- Due to this scaling of A the system-bath interaction is scaled with the 
dielectric constant e^. In the high temperature limit the reorganization energy is given by 

1113 

J{lo) 



\ = h 



duj- 



(4.43) 



/O 00 

This relation can be taken as motivation to scale the TB spectral density with eg like the 
reorganization energies A in the model with vibrational substructure as discussed in the 
subsection [4.2.4| . In the present calculations = = F is assumed. The absolute 



value of the damping rate F between the electronic states (see Fig. [4.5|) in this section 
is determined by fitting the FT rate for the solvent MTHF to be F = 2.9 x 10^^ s~^. It 
differs from the value given in section 4.5 because the different energies of states have 
been used here. 
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Figure 4.11: Variation of the potential minima for different solvents. Squares denote the 
bridge minima, circles the acceptor minima. The numbers correspond to the ordinal numbers 
in Table 4.3. The potentials are shown for solvent 6 (MTHF). 



The advantage of the TB model in comparison to the model with vibrational substruc- 
ture is the possibility to determine the transfer rate k^T and the final population of the 
acceptor state either numerically or analytically. For all situations described in this paper 
the differences between analytic and numerical results without the extra assumptions are 
negligible as shown in section [4.5|. 



In Fig. |4.11| it is shown how the minima of the potential curves change with varying the 
solvent due to the changes in Eqs. ( |4.37|) to ( |4.39| ). The solvents are listed in Table ^ 



together with their parameters and the results for the FT rates in both models. For 
larger eg the coordinates of the potential minima of bridge and acceptor increase while 
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Table 4.3: Parameters and obtained transfer rates for different solvents. The references behind 
the names of the solvents cite the sources of eg and eoo- T denotes the damping rate in the TB 
model. The ET rate for the solvent MTHF has been used to fix the damping parameter of the 
models. The reaction rates k^r^ were obtained using Eq. ( 4.32| ) within the TB model and the 
reaction rates within the vibronic model. 



solvent 












es 


f oo 


AG21 


AG31 


^21 


■^31 


r 


Ucl 
ft-ET 


iLvib 




















[eV] 


[eV] 


[eV] 


[eV] 


[10" s-i] 


[10^ 


[108 s-1] 


1. 


CYCLO 1140 








2.02 


2.00 


0.976 


0.393 


0.007 


0.012 


0.042 


0.181 


0.7 


2. 


toluene 


187| 








2.38 


2.24 


0.867 


0.202 


0.039 


0.069 


0.227 


1.04 


0.8 


3. 


anisole 


188 










4.33 


2.29 


0.590 


-0.281 


0.300 


0.524 


1.751 


4.24 


2.30 


4. 


dibromoethane 


188] 


4.78 


2.37 


0.558 


-0.336 


0.312 


0.544 


1.817 


4.63 


2.45 


5. 


chlorobenzene |187 


1 


5.29 


1.93 


0.529 


-0.388 


0.481 


0.839 


2.804 


3.21 


3.63 


6. 


MTHF |140| 








6.24 


2.00 


0.486 


-0.462 


0.497 


0.868 


2.900 


3.59 


3.58 


7. 


methyl acetate | 


18 




6.68 


1.85 


0.471 


-0.489 


0.571 


0.996 


3.328 


2.96 


4.15 


8. 


trichloroethane 


18 




7.25 


2.06 


0.454 


-0.512 


0.508 


0.887 


2.960 


3.98 


3.50 


9. 


dicliloro methane 14(J|| 


9.08 


2.03 


0.413 


-0.590 


0.559 


0.977 


3.264 


4.00 


3.80 



their energies decrease with respect to the energy of the donor. The energy difference 
between donor and bridge decreases with increasing eg. This makes a charge transfer 
more probable. For small eg the acceptor state is higher in energy than the donor state; 
nevertheless there is a small ET rate due to coherent mixing. For fixed e^o the ET rate is 
plotted as a function of the dielectric constant eg in Fig. |4.12| . The ET rate in the vibronic 
model increases strongly for small values of eg while the increase is very small for eg in 
the range between 5 and 8. The increase for small values of eg is due to the fact that 
with increasing eg the minimum of the acceptor potential moves from the position higher 
than the minimum of the donor level to the position lower than the donor level. So the 
transfer becomes energy favorable. This can also be seen when looking at the results for 
the TB model without scaling the electronic coupling with the Franck-Condon factor. In 
this case the ET rate increases almost linearly with increasing eg. The effect missing in 
this model is the overlap between the vibrational states. If one corrects the electronic 
coupling in the TB model by the Franck-Condon factor of the vibrational ground states 
as described in Eq. ( [4.42| ), good agreement is observed between the vibronic and the TB 
model. 

The ET rate for the vibronic model shows some oscillations as a function of eg. This is 
due to the small density of vibrational levels in this model with one reaction coordinate. 
All three electronic potential curves are harmonic and have the same frequency. So there 
are small maxima in the rate when two vibrational levels are in resonance and minima 
when they are far off resonance. Models with more reaction coordinates do not have this 
problem nor does the simple TB model. If these artificial oscillations would be absent, the 
agreement between the results for the TB and the vibronic model would be even better, 
because the rate for the vibronic model happens to have a maximum just at the reference 
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Figure 4.12: Transfer rate as a function of the dielectric constant Cs for both models together 
with experimental results |140| ]. The rates for the vibronic model are given by the circles. The 
dashed line shows the rate for the TB model with electronic couplings Vmn as in the vibronic 
model. The solid line represents the rate for the TB model with Vmn scaled as given in Eq. 



point eg = 6.24 which we have chosen to fix the spectral density, i. e. for MTHF. 

The comparison of the two models has been made assuming that the scaling of energies 
as a function of the dielectric function is correct in the Marcus theory. There have been 
a lot of changes to Marcus theory proposed in the last years. Marcus theory assumes 
excess charges within cavities surrounded by a polarizable medium and there one only 
takes the leading order into account. Higher order terms are included in the so called 
reaction field theory (see for example ||175|| ). But to compare different solvation models 
is out of the range of the present investigation. Some more details on this issue for the 
TB model are already given in section T^. Here we just want to note that the effect of 
scaling the system-bath interaction with eg, as assumed in the present work for the TB 
model, has no big effect on the ET rates. 



4.8 Summary 

We have performed a study of the ET in the supermolecular complex H2P — ZnP — Q 
within the TB model treated with the DM formalism. The determined analytical and 
numerical transfer rates are in a reasonable agreement with the experimental data. The 
superexchange mechanism of ET dominates over the sequential one. We have investigated 
the stability of the model varying one parameter at a time. The qualitative character of 
the transfer is stable with respect to a local change of system parameters. It is determined 
that the change of the dominating transfer mechanisms can be induced by lowering the 
bridge state energy. The physical reasons of system parameters scaling as well as the 
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relation of the theory presented here to other theoretical approaches to ET have been 
discussed. 

The validity of the TB model and its applicability to H2P — ZnP — Q is confirmed be- 
cause one gets good agreement for the ET rates of the models with and without vibrational 
substructure, i. e. the vibronic and the TB model, if one scales the electronic couphng 
with the Franck-Condon overlap matrix elements between the vibrational ground states. 
The advantage of the model with electronic relaxation only is the possibility to derive 
analytic expressions for the ET rate and the final population of the acceptor state. But 
of course for a more realistic description of the ET process in such complicated systems 
as discussed here, more than one reaction coordinate should be taken into account. 

The calculations performed in the framework of the present formalism can be extended 
in the following directions: (i) Considerations beyond the kinetic limit. The solvent 
dynamics has to be included into the model as well as, probably, non-Markovian RDM 
equations, (ii) Enlargement of the number of molecular blocks in the complex, (iii) Initial 
excitation of states with rather high energy should open additional transfer channels. 



Chapter 5 
Mixture of Solvents 



In this chapter we apply the RDM method which was described in chapter ^ to another 
experimental system i. e. porphyrin triad complexes. This method provides a quantitative 
description of time-resolved and steady-state properties such as fluorescence quenching of 
the porphyrins. It will be shown that our calculations do agree with already performed 
experiments. 

We start this chapter with a brief review of an important role of the supramolecular 
porphyrin array in electronic energy transfer, charge transfer, and photoinduced struc- 
tural changes in biological systems since such processes constitute the basis for the design 
of molecular devices of applicative interest and for the understanding of the most impor- 
tant photobiological processes. First, in section |5.1| we introduce the chemical structure 
of the self-assembled porphyrin triad ZnPD — H2P. Next, in section |5]^ we qualitatively 
describe the photophysics of the porphyrin triad in different solvents with various di- 
electric constants. In section |5.3| the parameters of the RDM formalism are discussed 
and determined for both the porphyrin complexes and the quantitative description of the 
fluorescence quenching. The explanation of the fluorescence quenching is done in sub- 
section |5.4.1| . The analysis of our calculations and the experimental results performed 
in subsection |5.4.2| yield valuable information about the reaction mechanisms. Finally, 
in section ^.5| we summarize the achievements of this chapter. Some generalizations are 
transferred to the Appendix O and require more detailed consideration in the future. 



5.1 Self-assembled triad of porphyrins 



Because of their widespread occurrence in photosynthetic reaction centres and other 
electron-transfer systems supramolecular porphyrin arrays have played a leading role in 
the study of energy and charge transfer processes in biological systems [|18 



19C, 191, 192 



1931 , |194|] or to gain insight into the principal possibilities of molecular electronics ||1 95 1 



Among supramolecular porphyrin arrays three-component covalently-linked donor-acceptor 
systems ||196| , |197|| have attracted a lot of interest, especially, in the context of photochem- 
ical molecular devices ||198|, |199||. Besides, it is possible to form donor-acceptor systems 



using both covalent bond and non-covalent binding self-assembling based on coordination 
interactions of metallo-porphyrins with appropriate extra ligands pUO| , |20 1| , |202| . These 
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systems are formed from chemical dimer of Zn-octaethylporphyrin with a phynyl spacer 
(ZnPD) and dipyridyl-dipheny [porphyrin with nitrogens in meta-positions of pyridil rings 
(H2P) via two- fold coordination of two central Zn ions of the dimer with pyridil rings of 
the extra-ligand. The chemical structure of a self-assembled porphyrin triad is represented 
in Fig. pTT] . It was experimentally found in these systems | |201|| that the population of 




Figure 5.1: Schematic presentation of ZnPD — H2P. 



the excited state of H2P is dependent on the solvent dielectric constant as well as on 
temperature. This is attributed to a close lying charge separated state which exchanges 
rapidly with the exited state of H2P. 



5.2 Photophysics of a porphyrin triad 



As determined by fluorescence spectroscopic measurements [p01|] , where one excites ZnPD, 
the triadic aggregate in non-polar solvents (toluene, methylcyclohexane) demonstrates flu- 
orescence quenching for ZnPD and quite intensive fluorescence for H2P. It is reasonable 
to assume that energy transfer processes may be involved in this quenching. The fluores- 
cence ZnPD — H2P* —>■ ZnPD — H2P (characterized by the fluorescence time rp = 7.7 ns) 
is less intensive with respect to that for pure H2P* H2P (tf = 9.3 ns) at the same 



conditions ||201|| . The quenching of H2P fluorescence increases with the polarity of the 
solvent. To check whether this is the case an experiment has been performed [p01|] . In this 
experiment the triadic complexes were formed in a solution of pure toluene (low dielectric 
constant) and then in an admixture of acetone (up to 20%) added to toluene (relatively 
high dielectric constant). In the first case the aggregate shows fluorescence with the mean 
band maximum of 716 nm which is attributed to H2P and fluorescence excitation spectra 
clearly demonstrate the presence of an energy transfer. While in the second case H2P 
shows noticable fluorescence quenching. Thus the excitation is lost without radiation and 
ET must take place. 
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Figure 5.2: Scheme of the photoinduced processes. 



The physics of this reaction is quahtatively shown in Fig. |5.2| , where D and A represent 
the electron donor ZnPD and the electron acceptor H2P, respectively. Note that here D 
and A denote other compounds than in chapter 4 while the symbols of plus, minus, and 
star are defined in the same way as in section 4.2. Light excitation of ZnPD to its singlet 
state |1) is followed by weak dissipative energy transfer to the exited state |3) (process a). 
Inaddition to this process, an ET (process b) takes place, leading to the charge-separated 
state |2). The evolution of the population of this state |2) is distributed between charge 
recombination (does not enter into the scheme) and further weak coherent and dissipative 
charge transfer (process c) to give the state |3). Because the decay of the charge-separated 
state occurs quite slowly for porphyrin- type systems ||203|| , we do not consider this process 
and do not include it in Fig. The energy of state |2) depends on the dielectric constant 
of the solvent. 



5.3 Density matrix model parameters 

The RDM model is able to describe coherent dynamics, dissipative dynamics, and ther- 
mally activated processes, thus the RDM formalism can be used for a quantitative ex- 
planation of the reaction which was described in subsection 5.2. We apply the system of 
equations (f4.20|) - (|4.21|) to describe the processes in the porphyrin triad. 



To calculate the population of the excited states physically reasonable values of the 
model parameters, such as coherent and dissipative couplings and energy of states have 
been chosen. The physical behavior of the system is determined in leading order by the 
energies, remains rather stable in respect to change of other parameters. Nevertheless 
it is reasonable to give some physical foundation in which region the parameters values 
should lie. 

In accordance with the consideration given in section 4.7 the parameters for the model 
with only electronic states can be extracted using some properties of the appropriate 
model with vibrational substructure. In such model one describes the porphyrin triad 
with the relevant potential energy surfaces in the space of a single reaction coordinate 
which reflect the extent of the solvent polarization induced by the field of the triad. In 
order to allow the analysis of absorption and emission spectra one should add to the model 
a ground state potential |0). For the neutral excited states |1) and |3) the difference in 
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absorption and emission spectra of the transition ground-excited determines the shift qio 
of the equilibrium point of the excited state \i) in respect to the equihbrium point of the 
ground state |0). The shift of the equihbrium point of the state |2) with charge separation 
can be calculated using the reorganization energy Eq. (4.39). 

We assume that the sign of the shift for the state |3) should be negative in respect 
to the signs of the shifts of |1) and |2). To check whether it is really so one should 
analyse the dependence of |1) |3) transition rate on the solvent polarity. Suppos- 
ing that the vibrational excitations again do not play a role and in accordance with 
the explanation after the Eq. (4.42) the coherent couplings can be expressed as Vi2 — 
Vi2exp[-{2hLJyib)-^{q2o - qiof], V23 = V23 exp [-(2/iu;vib)"^(g2o - ^so)^] etc.. Therefore 
the relation of the coherent couphngs reads 

^^12/^^23 = Vi2/V^23exp |-(2;ia;vib)"Mgio - qIo - 2g2o(gio - ?3o)]} ■ 

For qio ~ —qso and 520 3> —^30 the exponential reads exp [— (2/ia;vib)^^ 4g2o9io]- This is 
an exponential of a positive argument so the V12 should be at least e times larger as V23. 



Table 5.1: Coherent and dissipative couplings between the electronic states. 



Coupling Value, meV Physical Process Comment 



V12 
V32 

Vl3 



r 



12 



r32 



13 



60 electron transfer 

D*A D+A- 
3 hole transfer 

DA* D+A- 

12 energy transfer 

D*A DA* 

0.41 loss of co- 
herence for 
D*A D+A- 

2.50 loss of co- 
herence for 
DA* D+A- 

0.37 loss of coherence 
for D*A DA* 



induced by the wavefunction overlap 

weakened by the screening field of the 
electron from the LUMO of the accep- 
tor 

induced by the dipole-dipole inter- 
action of the excited states ^13 ~ 
Pd*aPda*/'^da^ 

interaction of the transition dipole mo- 
ment with environmental dipoles 

induced by the interaction with the en- 
vironment 



estimated by taking into account other 

dissipation 'jij = J2{dik + dkj) 

k 



The electronic matrix elements V^i, should also have a difference because the |3) — > |2) 
process imply the ET from one inner orbital HOMOd to another inner orbital HOMO a • 
So for the first attempt we take vu = 0.06 eV as typical value of the coupling be- 
tween porphyrins (same as the coupling (H2P* — ZnP — Q|if|F2P"'' — ZnP~ — Q) for zinc- 
porphyrin-quinone complex in chapter 4 with the precision up to the first significant figure) 
and V32 — 0.003 eV is taken to be 20 times smaller. 
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The range of reasonable values of relaxation constants T^^, can be estimated using 
Eq. (4.43) and reorganisanization energies from the model with vibrations. The values of 
the couplings are represented in Tab. 



5.1 



The precise determination of the parameters 
needs a special investigation within quantum chemical calculations. This is the next step 
in the theoretical description of this porphyrin triad. The quantum chemical calculation 
of the couplings should be done in the future. 
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Figure 5.3: Energy difference between states |2) and |3) 

The excited state energies -E'd*a = 2.10 eV and -Eda* = 1-91 eV are found spectroscop- 
The energy of the state |2) needs more attention. It depends on the solvent 



ically 1 140 



and can be calculated in any solvent with the help of Weller's formula [ ^04|| : 



D+A- 



Eb+a- (ft. 



1 

et 



+ 



1 



Aireo V2rD 2rA Tda 



(5.1) 



where e denotes the solvent static dielectric constant, td = ta = 5.5 A are the donor 
and acceptor radius, respectively, tda = 8.8 A is the distance between them. In our case 
the solvent consists of the main compound toluene with a dielectric constant et and a 
small concentration c of the additional compound acetone with a dielectric constant e^- 
In accordance with 



205 the effective dielectric constant of the mixture reads 



3(e« 



et + c- 



etjet 



ea + 2et 



(5.2) 



The change of e induces an energy shift of the state |2). We present the energy difference 
AE = E^ — E2 between states |3) and |2) as a function of c in Fig. [5.3|. The values 
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Eb+a- = 1-90 eV, ct = 2.38 for the pure toluene solution and ea = 10 for the acetone were 



defined in 201 



5.4 Physical processes 

5.4.1 Fluorescence quenching: simulations and experiment 

In the porphyrin triad the competition between charge transfer and energy transfer (pro- 



cess a and b in Fig. |5.2| , respectively) cause a rather complex dynamics. Mathematically 
it is easy to calculate the system state at the infinite time t = oo. For our real system the 
electron transfer time tet is much shorter than the time of fluorescence rp. Thus we ap- 
proximate t = oo with some time moment when the ET has finished and the fluorescence 
has not occured yet. On this time tet < t < Tp the system reaches the quasi thermal 
equilibrium between the excited state |3) and the close lying charge separated state |2). 
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Figure 5.4: Solvent induced fluorescence quenching, T 
rescence quenching, Ca = 7% (right). 



293 K (left); freezing induced fluo- 



We have calculated the equilibrium population of the state |3) numerically with the 
RDM-method. This population corresponds to the presence of the fluorescence into the 
ground state. We denote the population of this state with Pda* = Psa- 

It has been found that ^33(00) decreases in two cases: (i) lowering of the energy i^D+A- 
induced by increase of aceton concentration, see Fig. |5]^ (left) and (ii) lowering of the 
temperature, see Fig. |5.4| (right) This decrease corresponds to the experimentally observed 
fluorescence quenching P01|| . 



5.4.2 Reaction mechanisms 



The variation of the acetone admixture concentration and temperature changes the flu- 
orescence intensity P33 as well as a character of the excited states dynamics. Our sim- 
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ulations display that the time dependence of the population pss (Fig. |5]^) as well as its 
temperature dependence (Fig. |5l^ ) at low acetone concentration quahtatively differs from 
the behavior at high acetone concentration. This is the most important result of the 
RDM-calculations. It demonstrates the qualitatively different reaction mechanisms under 
various experimental conditions. 




Figure 5.5: Influence of concentration on the reaction mechanism, T = 293 K, at low Ca 
the reaction passes in one step (circles and diamonds), high Ca induces the two-step reaction 
mechanism (triangles, squares, and stars). 

A low concentration of acetone induces a low energy detuning between the states |2) 
and |3). In this case (see Fig. 5.6 for smah Ca ), pss starts to increase with time due 
to energy transfer (process a) and then does not change as |3) reaches the quasi-thermal 
equilibrium. Thus the equilibrium population p33(oo) is reached in one-step and a reaction 
rate k can be found with an one-exponential fit Eq. ( [4.28| ). 

As shown in Fig. in the case of high aceton concentration the energy detuning 
between states |2) and |3) becomes larger and in addition to the energy transfer (process 
a), the hole transfer (process c) takes place, thus the equilibrium population pss is reached 
in two steps. At first, the energy transfer creates a time-dependent maximum of the 
population of |3) (see Fig. 5.6 for large Ca) and then hole transfer |3) |2) slowly induces 
depopulation of |3) down to the equilibrium population. The reaction occurs with the 
help of a sequential transfer, which is described by two rates (increase and decrease) and 
the one-exponential fit for k cannot be used in this case. 

The temperature dependence of p33(oo) (Fig. 5.7) also refiects the change of the re- 



5.5. SUMMARY 



71 




Figure 5.6: Thermally activated mechanism, the population remains constant until energy 
difference is lower than the thermal energy. 

action mechanism which depends on the acetone concentration and temperature. For 
a low acetone concentration the equilibrium population ^33(00) increases with tempera- 
ture because AE between states |3) and |2) is lower than the thermal energy. This case 
corresponds to a one-step reaction. 

In the case of a high acetone concentration AE becomes larger. Thus we have two 
regimes ksT > AE and ksT < AE. The first case correspond to the one-step reaction, 
the second to the two-step reaction. For high acetone concentration the increase of tem- 
perature induces the crossover from the second to the first type of behavior: While the 
temperature is quite low, its increase does not change ^33(0x3). It occurs because in the 
abscence of the thermally induced transitions ^33(00) is determined in leading order by 
the coherent mixing V23 and does not depend on the temperature. Then, when the tem- 
perature is quite high in comparison to AE, so that ksT ~ AE, the reaction crossover 
into the one-step regime is detected and P33 starts to grow with the temperature. For a 
high acetone concentration (20%) there is no essential temperature dependence of P33 in 
the considered interval of temperatures. 



5.5 Summary 

The RDM method explains the physical properties of a real system, e. g., a self-assembled 
porphyrin triadic aggregate both qualitatively and quantitatively and gives good agree- 
ment with already performed experiments. However, in order to describe the experiments 
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completely it is necessary to take into account a pump-probe laser field. This general- 
ization is quite complicate and requires more detailed consideration in the future. We 
represent some mathematics involved in this generalization in Appendix |^. 

Our calculations predict a solvent variation induced crossover between two regimes of 
the time evolution that most probably follows from the energy level dependence. There 
is no experimental data about this phenomena, that is why it is interesting to prove the 
existence of the described crossover in experiments. 



73 



Chapter 6 
Conclusions 



It is shown that the vibrational wave packet relaxation of initially coherent (displaced) 
states as well as the quantum superposition of coherent states in heat baths with different 
spectral densities exhibit a number of peculiarities compared with the cases of linear and 
quadratic systcm-bath interactions. A strong dependence of the relaxation rate on the po- 
sition of the spectral density maximum of the bath is found. The difference discriminates 
the mechanisms of the molecule-environment interaction. 

Based on the RDM method, we calculate the dynamics of ET for systems consisting of 
donor, bridge and acceptor in different solvents. In our first approach it is assumed that 
vibrational relaxation is much faster than the ET. Transfer rates and final populations of 
the acceptor state are calculated numerically and in an approximate fashion analytically. 
In wide parameter ranges these solutions are in very good agreement. The theory is 
applied to the ET in H2P — ZnP — Q with free-base porphyrin (H2P) being the donor, 
zinc porphyrin (ZnP) the bridge, and quinone (Q) the acceptor. It is shown that the 
transfer rate can be controlled efficiently by changing the energy of the bridge level that 
can be done by changing the solvent. The effect of the solvent is determined for the 
models of single and multiple cavities. This approach has been compared to the second 
approach, where a vibrational substructure is taken into account for each electronic state 
and the corresponding states are displaced along a common reaction coordinate. In both 
approaches the system is coupled to the bath of HOs but the way of relaxation is quite 
different. For the comparison of these two models the parameters are chosen as similar 
as possible for both approaches and the quality of the agreement of the approaches is 
discussed. 

Applying the RDM theory to the photoinduced processes in ZnPD — H2P it has been 
found that the population of the state ZnPD — H2P* which controls the intensity of fluo- 
rescence of this complex is strongly influenced by the temperature and dielectric constant 
(polarity) of the solvent. The change of the last two parameters alters the character of 
the dynamics of the state ZnPD — H2P*. 
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Chapter 7 
Outlook 



The work is devoted to the investigation of the influence of a heat bath on the physical 
processes in a quantum system. We use the density matrix theory as one of the most 
powerfool tool for investigation of quantum relaxation. In the beginning of the work we 
mention and recall the most important steps of derivation of the equation of motion for 
the RDM (master equation) for an arbitrary quantum system in diabatic representation 
interacting with the environment modeled by a set of independent HOs. 

At first we apply the theory to a single state in the diabatic representation decoupled 
from other states but having vibrational substructure presented by a single HO. We have 
performed a thorough investigation of this model with the help of the master equation, 
which has been solved analytically and numerically. For this system the wave packet 
dynamics in coordinate representation has been analysed for two models of the bath and 
two initial states. The different models of the bath have their maxima of the spectral 
density near the system frequency and near the double of the system frequency. The 
considered initial states are a coherent state and a superposition of coherent states. It has 
been shown that the wave packet dynamics demonstrates either " classical squeezing" and 
the decrease of the effective vibrational oscillator frequency due to the phase-dependent 
interaction with the bath, or a time-dependent relaxation rate, distinct for even and odd 
states, and partial conservation of quantum superposition due to the quadratic interaction 
with the bath. The decoherence also shows differences compared to the usual damping 
processes. There are two universal stages of relaxation which allows analytical solution: 
the coherence stage and the Markovian stage of relaxation. 

The density matrix theory for a quantum system in a diabatic representation has 
been applied twice to a study of the ET in the supermolecular complex H2P — ZnP — Q, 
namely with and without account for vibrations in the complex. With help of the model 
without vibrations we have determined analytical and numerical ET rates which are in a 
reasonable agreement with the experimental data. The superexchange mechanism of ET 
dominates over the sequential one. We have investigated the stability of the model varying 
one parameter at a time. The qualitative character of the transfer is stable with respect to 
a local change of system parameters. It is determined that the change of the dominating 
transfer mechanisms can be induced by lowering the bridge state energy. The physical 
reasons of system parameters scahng as well as the relation of the theory presented here 
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to other theoretical approaches to ET which do not accounts for the vibrations have 
been discussed. The vahdity of the model without vibrations and its applicability to 
H2P — ZnP — Q is confirmed because one gets good agreement for the ET rates of the 
models with and without vibrational substructure, if one scales the electronic coupling 
with the Franck-Condon overlap matrix elements between the vibrational ground states. 
The advantage of the model with electronic relaxation only is the possibility to derive 
analytic expressions for the ET rate and the final population of the acceptor state. 

We have also applied the RDM theory to explain the physical properties of a self- 
assembled non-fluorinated triadic porphyrin aggregate. We have simulated the processes 
in this complex placed in a mixture of solvents with different dielectric constants. Our 
simulations reproduce the intensity of the fluorescence of the aggregate and its dependence 
on temperature and mutual concentration of constituents in the mixture both qualitatively 
and quantitatively with reasonable agreement with already performed experiments. Our 
calculations also predict a solvent variation induced crossover between two regimes of 
time evolution. There is no experimental data about this phenomena, that is why it is 
interesting to proof the existence of the described crossover in experiments. 

The calculations performed in the framework of the present formalism can be extended 
in the following directions: 

1. Considerations beyond the kinetic limit. The solvent dynamics has to be included 
into the model as well as, probably, non-Markovian RDM equations. 

2. Enlargement of the number of molecular blocks in the complex. 

3. Initial excitation of states with rather high energy should open additional transfer 
channels. 

4. For a more realistic description of the ET process in such complicated systems as 
discussed here, more than one reaction coordinate should be taken into account. 

5. In order to describe the experiments completely it is necessary to take into account 
a pump-probe laser field. 

This generalization is quite complicated and requires a more detailed consideration 
in the future. In this work we have represented some mathematics involved in this 
generalization. 

6. Quantum chemical calculation of coupling matrix elements. 

The future development of the extension (1) lies in the correct treatment of non- 
Markovian effects in the description of the photoinduced charge- and exciton transfer of 
a single molecule in a solvent by the density matrix theory. The model should remain the 
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same as in this work although the solvent correlation functions should enter into the theory 
in explicite form without approximating them as a delta functions. For the calculation 
of the correlation functions of the bath of the polar molecules one should take the model 
of the dipole-dipole interacting spins, similar to the well known Heisenberg model. The 
response function of this spin lattice should be calculated for the time dependent charge 
separation in a single molecule embedded in this lattice. The proposed extension of this 
work would have the following advantages in respect with the traditional treatment of 
the bath: (a) It accounts for non-Markovian effects, (b) The state of the environment 
depends on the state of the system, (c) The solvent dipole value is the realistic parameter 
which enters into this model. As a possible outcome of the investigation of the bath of 
polar molecules modelled by the lattice of the rotators which interact as dipole-dipole one 
could expect some phase transitions with temperature as it often occurrs for interacting 
spins. 

As a general conclusion we could mention that coupling of a quantum system to the 
heat bath leads to the loss of energy, to disappearence of phase information, and ensures 
the irreversibility of the processes in the system. The influence of the environment should 

be in most cases included in the description of a real physical system. We believe that 
the developed methods and obtained results can be applied for other initial states and 
different couplings with the environment in real existing quantum systems. 
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Comparison with the 
Haken-Strobl-Reineker formahsm 



We compare the RDMEM ( |2.29[ ) with an analogous equation within the HSR model 

There is a term in the RDMEM of the HSR model, which is 



131, 132, 133, 178 



absent in our calculations. Here we show that this term is nothing but the difference 
between the full relaxation operator and the relaxation operator in RWA. We neglect this 
term corresponding to 7^,^ both in the equation of motion and in the expression for the 
transfer rate k^T due to the RWA. The symbol 7^1, is used in the HSR for the rate of 
changes in the system state induced by this term. 

First we mention that the Eq. (2.12) from ||178|| is derived under assumption 7^,^ = 0. 
The RWA relaxation term within our formalism Eq. ( p.29| ) yields the same RDMEM as 

I of this work without 7^jy. On 



Eq. (2.12) of Ref. ||178|| . This equation is solved in section | 
the other hand the transformation of non-RWA term into matrix form gives the expression 
associated with 7^,^ in the stochastic Liouville equations (SLE) formalism |]53| as we show 
here. 

The relaxation operator obtained within RWA Eq. ( ^.291 ) is assumed to describe the 
major contribution to the system dynamics. The importance of the non-RWA counterpart 



Lnon-RWAcr = ^ {r^,, [n{ujf,i,) + 1] ( V^u^, Vf,u 



+ 



+ 



+ 



)} 



(A.l) 



obtained using the interaction Hamiltonian Eq. (|4.14|) instead of Eq. (|4.16|) remains ques- 
tionable. 

The matrix form of Eq. ( |A.1| ) reads 



(-^non-RWA 0" ) k> 



{Eak [2n(t^AK) + 1] + r«A [2n(t^KA) + 1]} CTXk- 



(A.2) 



Thus, the full relaxation dynamics for the system-bath coupling Eq. ( [4.14|) is described as 
the sum of two terms: Eq. ( [4.18| ) and Eq. ( [A.2| ). To compare the present approach with 



the SLE formalism used by HSR ||131| , |132|| we recall the application of both methods to 
the simplest system, namely the TLS. The SLE method |133[ provides the following 
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structure of the incoherent term for the RDMEM (Eqs. (3.1a)-(3.1d) in |53[]): 
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Applying the present model to the TLS we obtain 



(A.3) 
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(A.4) 
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(A.5) 



The relaxation coefficients 71 and 7_i from Eq. ([A .31) that are, evidently, the difference, 
mix the non-diagonal elements of DM (T21 and au. 

The derivation within the formahsm of present paper without RWA Lrwa + -^non-RWA 
ensures the same structure of relaxational dynamics as Eq. ( |A.3D . In details: the non- 
RWA terms r[2n(ti;2i) + 1] from Eq. ( |A.5| ) ensures the same effect as coefficients 71 and 
7_i from Eq. ( |A.3D . But such kind of terms is really absent in the RWA relaxation term 
Eq. ( [A.4|) which we use for the material-oriented calculations. 

At the present stage it is possible to estimate that consideration of non-RWA terms 
makes a smooth change of the characteristics of the process, i.e., /cet while the expressions 
lose their simplicity. In our opinion the desired comparison of the more (without RWA) 
or less (with RWA) precise description of the TB model goes beyond the goals of this 
work. 

Another difference of equations (|A.3|) and ( |A.4|) is that the first of them (SLE method) 
describes the same dissipative transition probability from (T22 to an and back. It leads 
finally to the equal population of both levels. In our approach it is possible only if the 
states |1) and |2) are isoenergetic. In any other case the transition from upper level to 
the lower one will be more intensive as an inverse transition to construct the Boltzmann 
distribution of populations at the infinite time. 

The third, and, perhaps, the main difference of SLE and our methods is that SLE 
method assumes the modulation of system frequency 



(A.6) 



(Eq. (2.10) in Ref. H], Eq. (2.5) in Ref. [|13T|], Eq. (2.2) in Ref. [|3 



where de- 



notes exciton annihilation operator, h^y{t) stochastic function {h^y{t)) = 0, while in our 
approach it is assumed that system performs exchanges of quanta with the quantized 
modes of the thermal bath. So we conclude that the used RDMEM within TB model 
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coincides with well-known HSR equation for exciton motion under certain approxima- 
tions: i) energy levels are isoenergetic for the first RDMEM and ii) RWA for the second 
one. The similarity of the equations appears although different models are used for the 
environment. The generalization ||206|| of the SLE method appeals to the quantum bath 
model with SB coupling of the form H^^ 
transition frequency. In Ref. 



206 



o-x + ctxj , which modulates the system 
the equations for exciton motion are using projection 
operator technique without RWA leading to the presence of 7 in the generalized stochas- 
tic Liouville equation (GSLE). So, taking the different SB coupling we have rederived a 
RDMEM which coincides with the GSLE |20|1 after applying RWA. Both GSLE and our 
RDMEM are able to describe finite temperatures and non-periodic systems. 
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Appendix B 

Full Model: Vibrations, Optics, 
Memory Effects 



In the full model the molecule is irradiated by the electromagnetic field as sketched in 
Fig. p.l| . The common Hamiltonian ( |2.1| ) is extended by the field term H^^ and is written 




Figure B.l: Schematic presentation of full model, 
as follows H = + + H^^ + H^^. Here the diabatic system Hamiltonian 



H'' = hY,HfjI,n){I,n\+ Vi,n,J,m\I,n){J, 

I I,n,J,Tn 



m\ 



is characterised by the electronic- vibrational matrix elements iff ^ = Qj + uji{n + 1/2) 
and coherent mixing matrix elements V/,n,j,m = VuF^ci,! .''^^ J, m)(l — 5ij). Here VLi and 
uji stand for electronic and vibration transition frequencies, Vu for electronic coupling, 
Fpc{I,n, J,m) = {I,n\J,m) for Franck-Condon factors. The system couples to the bath 
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of HOs as follows H^^ = [E Ki{cj + ci)] [E k^ibt + b^)], where c/ = E Vn\I,n-l){I, n\ 

I 5 n 

stands for annihilation of a vibronic quantum in the Jth electronic state. The electro- 
magnetic field E{t) = Eq e^^^ + Eq e~^^^ is described by frequency W and strength Eq = 
Eq = Eq. The system-field interaction is written as H^^ = E \Eo\pi^rn,j,n\I , n) {J, m\ 

I,n,J,m 

with the dipole pi,n,j,m = e(/,n|E + ci)\J,m). Taking the notation W/,n,j,m = 

Wsign{Hf^^ - Hl^) we obtain the relevant RDMEM: 

,171, K,l 

+ Vi 

,171, K, l){PK,l,I — Pl,m,K,l) + Li 

,m,I,m 

K I 

Pl,m,J,n = [-i{Hf „^ - Hj J + iWl^rn,J,n]pi,m,J,n (B.l) 

,m,J,n + Vj^rn,J,n){Pj,n,J,n — Pi ,m,I,m) ~r J-'I ,m,,J,n 

Li,m,j,n= - (m[7,, + 23f?7a + 2^^^j] + 23f?7yj)p,,^,j,^ 

,J,n+l 

,m— 1, J,n— 1 • 

The relaxation functions 77J are derived from the bath correlation functions: 

Kik^Kjk^N^LJ^) exp {—iu^t + iu^r + iujjt — iuijr) , 
in the same way as introduced in chapter |^. 
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Thesen zur Dissertation 



• Die analytische Losung der Bewegungsgleichung fiir die reduzierte Dichtematrix 
wurde fiir einen harmonisches Oszillator und fiir die Kopplung an ein thermisches 
Bad fiir Zeiten kiirzer als die Badkoharenzzeit und fiir lange Zeiten berechnet. 

• Die gefundene Wellenpaketdynamik in Koordinatendarstellung fiir ein Bad mit Spek- 
traldichtemaximum bei der Systemfrequenz zeigt eine Oszillation des Koordinaten- 
mittelwerts und eine oszillierende Verringerung der Breite des Wellenpakets und der 
Energie des anfanglich koharcnten Zustands, wobei die Koharenz der Superposition 
von koharente Anfangszustanden schnell abklingt. 

• Die Simulation fiir Modelle mit Badfrequenz, die doppelt so hoch wie die Systemfre- 
quenz ist, zeigt eine zeitabhangige Rate des Energieverlustes, die sich fiir gerade und 
ungerade Oszillatorzustande unterscheidet, mit einer relativ langen Dekoharenzzeit 
eines Superpositionszustandes. 

• Fiir die Losung einer Dichtematrixbewegungsgleichung in diabatischer Darstcllung 
und unter Vernachlassigung der vibronischen Struktur fiir angeregte Zustandc des 
H2P— ZnP— Q -Komplexes wurde gefunden, daB die Besetzung des ladungsgetren- 
nten Zustandes H2P''" — ZnP — exponentiell bis zum Gleichgewichtsverteilungs- 
wert wachst. 

• Der Superaustausch-Transfermechanismus ist verglichen mit dem sequentiellen Trans- 
fer im Komplex H2P— ZnP— Q vorherrschend. Dieser SchluBfolgerung bleibt auch 
bei lokaler Anderung eines Systemparameters giiltig. 

• Die Verringerung der H2P^ — ZnP~ — Q — Zustandsenergie ruft die Anderung des 
Transfermechanismus von Superaustausch zu sequentiellem Transfer hervor. 

• Bei der Behandlung der dielektrischen Umgebung in einem Modell, das entweder 
Donator, Akzeptor und Briicke gemeinsam im einem Hohlraum des Dielektrikums 
enthalt oder jeden Baustein in einen einzelnen Hohlraum, zeigte sich, daB die Trans- 
ferraten im zweiten Fall genauer beschrieben wurden. 

• Die Simulation des Elektronstranfers unter Beriicksichtigung der Schwingungen ergibt 
Transferraten, die mit den analytischen Transferraten vom Modell ohne vibronische 
Unterstruktur gut iibereinstimmen, woraus wir schlieBen, daB schwingungslose 
Modelle geignet sind, Elektronentransfer in H2P — ZnP — Q zu beschreiben. 
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• Die Fluoreszenzstarke des selbst-aggregierten molekularen Komplex ZnPD — H2P 
in einer Mischung von Losungsmittcln mit vcrschicdcncn diclcktrischcn Konstanten 
wachst mit der Temper atur und sinkt mit steigender Azetonkonzentration. 

• Die Simulationen sagen einen von Losungsmittel induzierten Ubergang des Reak- 
tionsmechanismus im Komplex ZnPD — H2P voraus. 
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